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Some Physical Properties of Porcelains in the Systems 
Magnesia-Beryllia-Zirconia and Magnesia-Beryllia- 
Thoria and Their Phase Relations’ 


By S. M. Lang, L. H. Maxwell, and R. F. Geller 
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ternary eutectic for the system MgQO-BeO-ZrO, was located at mole composition 
IMgO:5BeO:3ZrO,y and 1,669° C, and that for the svstem MgQ0-BeO-ThO, was located at 
mole Composition 5MeQ 1OBeO: lL ThO, and 1.707° ¢ No binary or ternary, compounds 
were found in either of the systems, and an area of solid solution is indicated for eacl \ 


number of bodies were found suitable for various high-strength and high-temperature services 


where thermal shock resistance is not of prime importance The compressive strength 
values at room temperature ranged as high as 266,000 Ib in? for some porcelains in both 
system The maximum values for the strength in bending at room temperatures and at 
1Sso0° | WRS2° ¢ were considerably higher for the MgO-BeO-ZrOQ. bodies (39,000) and 
$2,000 Ib in, respectively) than for the MgO-BeO-ThO, bodies (14,000 and 16,000 |b in 


respectively), and the modulus of rupture va 


binary system MgQO-ZrQ. 20 0000 Ib in 


Ve0-BeO-ZrO 
All of the 


Was 


porcelains was fair to good, a 
poor 
inn thre 


porcelain industry, whereas part of the 


MgQ-ZrQ, porcelains can be so matured 


I. Introduction 


Under the impetus of the recent war and because 
of the outstanding advances made in the design 
and efficiency of such power plants as the gas- 
turbine and jet-propulsion engines, varied pro- 
grams are being sponsored by industrial, research, 
and government laboratories in the realm of high- 
temperature materials 

In the 


present models of the various heat power-plants 


design and, especially, in’ efficiency 
are severely restricted by the ability of some of 
their components to withstand the required high 
The im- 


portant physical properties of such power-plant 


operational temperatures and stresses 


elements mre dependent to on large extent upon 


Whether the part in question ts to be stationary 


or moving, Whether the eagine is to be continu- 
ously ol intermittently operated, or any com 
bination ol these conditions iI] 

on N ki 


NA 7, 1 ‘ N NKR 


Refractory Oxide Porcelains 


MeO-BeO-ZrO. bodies can be 


zw at 1.S00 i QRXx2° ( of one bodv in the 
The resistance to thermal shock for the 
nd for the MgO-BeO-ThO, porcelains it was 


matured in commercial kilns used currently 


MgQ-BeO-ThOs porcelains and none of the 


The allowable metal-blade temperature of the 
present conventional gas turbine, for example, is 


limited to approximately 1,500° F 


815° C). but 


conventional fuels could) provide inlet-gas tem- 


peratures exceeding 3,500° F (1.927° © 1] and, 
thereby, greatly increase the efficieney of such an 
engine Inaddit ion to refractoriness, considerable 


mechanical strength is required in certain in- 


Stances us illustrated, for example, by the tensile 
stresses set up ina turbine blade 
2, 3,4] 
required, 


by centrifugal 


force A number of 


other properties alse 
Which 


high 


resistance to oxidation and erosion 


are umong are: high thermal 


shoch resistance, thermal conductivity. 


and ability to 


Daintain exacting dimensions | high tem 


peratures 
Many of the ceramic oxides and combinations 
ol these oNxides re well suited lol such diversified 


uses as turbine nozzles, 


Interstave stator blades, 


rotors and buckets, rocket liners and nozzles. and 


eontre | Vahes because of there hioh melting points 


low density, and potentially low cost and availabil 
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ity Thi limitations imp mca yy the vlassy bond 
of the conventional feldspathic bodies and the 
earlier developments leading to the formulation of 


the glass-free oxic porcelains are discussed in 


detail inoa Excepting the 


cle tatled 


previous report 5 


studies of the phase relations of the 


oxides |6 » Many reports of the pure ceramic oxide 


and of the ceramu oxide-metal combination 
studies are still in the restricted classification 

In order to supply some of the fundamental 
information upon which the practical work of pro- 
ducing ceramics to meet the future needs can be 
based, there was undertaken a study of the charac- 
teristics ol variously heated binary and ternary 
combinations of the high melting-point oxides. A 
preliminary report was published in 1941 [7], and 
areport onsomeof the porcelain compositions of the 
systems MeOQ-BeQ-ALOs;, CaO-BeO-AlO,, BaQ- 
}6O, BeO-ALO,-ThO,, and BeO-AlLOs-ZrQ,, along 
with the probable phase relations in these systems, 
was published in 1946 [5]. A preliminary report 
on the strength in tension and the creep charac- 
teristics of a number of bodies described in the 
1946 [2], and the 
final report was published as NACA Technical 
Note L561, April IS, 1949. The present 
reports the results obtained since the 1946 report [5] 


1946 report also was made in 
papel 


lt ETI AS be of interest to the reader to compare 
some of the reported physical properties of a few 
bodies and of some. oxide 


refractory ceramic 


porcelain bodes. These data are given in table - 


II. Materials 


The oxides used in the preparation of specimens 
for testing the various physical properties were 
commercially available high-purity = materials, 
Whereas the oxides used in the study of the phase 
relations were of the highest purity available 

Beryllia (BeO): (1 


inal 99.7-percent purity 


Commercial grade of nom- 
The spectrograms show 
only traces of copper, iron, and magnesium and 


verv weak lines for silicon 


9) Fluorescent grade of over 99.9-percent 


purity 


Vagnesia (MgO 1) Fused, reervstallized ma- 


terial of nominal 97-percent purity 


2 Especially selected, water-clear crystals ol 


artificial periclase of over 99-percent nominal 


purity (furnished through the courtesy of The 


Norton Co 
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material of ove 
The 


show only very weak lines for calcium, copper, an: 


Calcined 
purity 


ThO l 


nominal 


Thoria 
99-percent spectrogran 
silicon 
2) Caleined material of at least 99 99-perce! 


purity 

Zirconia (ZrO Commercial grade of nomin: 
99-percent purity. After recaleining at 1,440 
C the spectrograms showed medium lines fo 


columbium, titanium, and only very weak lines o 
traces for calcium, copper, iron, magnesium, lea: 
and silicon. ‘This material was used for both th 


porcelain test) specimens and = for the phas 
relation studies 

As prepared for use, the materials were, in a 
Instances, sufficiently finely divided to piss t| 
No. 325 U.S. Standard Sieve In order to obta 
this desired fineness, it was necessary to redur 
the particle sizes of all of the materials except tl 
fluorescent grade bervilia. The commercial grad 
of bervllia was reduced to about 95 percent © 
minus $25 mesh material by wet grinding for 2 | 
ina porcelain mill with alumina balls (appros 
( 


mately 2% total mill gain of which about 1 
Both grad 


of magnesia could be reduced in size fairly eas 


was gained from the alumina balls 


ina mullite mortar, but it was found more ady 
able for quantity work to mill the less pure mat 
rial in aleohol for 30 min to 1 hr in a porcelain m 
with alumina balls (less than 's°% total mill gain 
Both the thoria and the zirconia were dry-grout 
in a steel rod-mill for 30 hr and then cleaned. A 
magnetic separator was used first, and the materi: 
was then washed alternately with 1:3 hydrochlos 
acid and distilled water until both the decant: 
acid and wash water gave a negative iron test usil 
potassium thiocyanate or potassium ferroevani 


or both as indicators 


III. Methods 


In order to insure the desired reproducibility 
results and to prevent segregation of the bat 
constituents during specimen preparation, a de! 
nite compounding method was followed. 

The calculations of weight composition from t! 
were based upon tt 


desired mole composition 


weights of the oxides ‘ 


0.005 weight perce! 


ae cepted molecular 


Calculations were made to 


All temperature values are given according to the Internatior I 


ture Scale of 192 
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and there was no correction factor applied for the 
percentage purity of the material at any time nor 
for anv ignition loss at 1.200° C when that loss 
was less than 0.01 ¢ of the various 30- to 400-¢ 
batches. Each oxide was weighed on a semi- 
analytical balance to 0.005 ¢ Approximately 
1.2 ml of distilled water per gram of batch was 
added to make a slip similar in consistency to a 
thick cream and one drop of a L0-percent water 
solution of Aerosol wetting agent was added for 
approximately each 10 ¢ of dry bateh The slip 
Wits machine blunged with a chemical laboratory- 
tvpe variable-speed mixer with plastic stirrers for 
not less than 15 min, damp-dried under infrared 
lamps, and completely dried in an oven at 110°C 
for at least 15 hr. Between 8S and 10 weight per- 
cent of a 5-pereent soluble starch solution was 
ndded as a binder and the dampened batch Wiis 
mixed and sieved «a number of times through a 
No. 30 [ S. Standard Steve lt was then pre- 
pressed at about three-fourths of the final forming 
presstuire to promote agglomeration, and the resul- 
tant briquets were crushed to pass either the No 
20 or the No. 30 sieve. The dampened mixture 
Was tamped and vibrated in the molds before the 
final pressing This pressing was done in case 
hardened steel, and in oil hardened, nondeforming 
tool steel molds liehtly lubricated with nu Wate 
emulsifving oil 

Specimens for the softening-range and fusion 
tests [5] were ground from by In. diameter 
plaques pressed at 15.000 Ib int. Those for the 


preliminary investigation of the maturing range 


ith high Iyy nh. Uhh diamete ! and for the cColn- 
pression tests (1'; in. high by in. in diameter 
were pressed at 9400 Ib in The thermal shock 
and stre neth test bars were pressed about nh 


thiek in a *y-in. by 6 in. mold at 6.500 Tb in 

The special high-temperature thoria-resistor 
furnaces used in these studies for determining 
Vitrification ranges, for maturing the compression 
specimens, and for the phase studies, have been 
deseribed (7, Ol, as have THhaenny of the test methods 
y These are re yp ated both for the convenience 
of the reader and to give the details of those 
procedures that have been revised 

driefly, the study of the phase relations was 
pee omplished by noting the melting characteristics 
of small, four-sided pyramidal specimens with an 
optical pvrometet and by petrographic and X-ray 
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examinations of the nonquenched heated spec 


mens. It must be kept in mind that equilibriun 
conditions at the elevated temperatures may no 
have been established [5] and that the phase rela 


tions suggested were based primarily on the fusion 


behavior observations. Along with the inherent 


difficulties of petrographic examination and = in 
terpretation of powder slides and thin-sections o 
these nonquenched specimens with high refractiv: 
indices (2.08 and 2.15 for ZrO, and ThO,, rv 


spectively), there were available only very smal! 


specimens for both the petrographic and X-ray 


examinations, 

Absorption determinations were made on spec 
mens boiled in carbon tetrachloride, and th 
results were converted to equivalent water absorp 
tion values. Mature bodies were those having 
less than one-tenth percent of equivalent wat 
absorption 

The length-to-diameter ratio of all compressior 
test specimens was approximately 2. All spec 
men ends were ground parallel on a Blanchare 
grinder using diamond wheels, and the specimer 


were crushed between parallel-ground, °.-in. thie! 


cold-rolled steel blocks ata loading rate of abou 
31,000 Ib in? per minute 

Following maturing of the thermal shock an 
rupture test bars in a high-temperature gas-fire: 
kiln, the specimens were ground on a Blanchars 
erinder using diamond wheels so that the bar 
would be approximately 'y-in. thick with parall 
faces. Ten eveles of quenching from 1,700° | 
931° () to an air-blast constituted the therma 
shock test. The arrangement of the equipme 
for this test is shown schematically in figure | 
“A. At! BY isan enlarged view of the specime 
holder with six test bars. The bars were plac: 
in the hot furnace (1,.700° F 5° BF), heated fe 
30 min (the furnace reached and maintains 
1.700° F for about the last 15 min), and we 
rapidly withdrawn to the air-blast. The oritir 
of the blower, delivering 625 ft ° of room-tempx 
ature air per minute, Was tin. in diemeter 2 
wis located on above the center of the vroup 
specimens. After remaining in the air-blast ff 
15 min (the specimens were only warm after abo 
5 min), each bar was examined for failure. Tl 
order and position of the test bars were thy 
changed according to a sehedule and were 


placed in the furnace to Devin the next evele 
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The modulus of rupture values in bending were 
both the 


determined at room temperature for 


quenched and nonquenched bars [10]. 
The modulus of rupture and modulus of elas- 
values at 1,800° F (982 


ticity (Young’s modulus) 











WEEE 
QM. 
SZ y 
m ¥ 


- 


nae 
ee ere aa 


























(") for the compositions studied in this tnvestiga- 
tion were determined by using a furnace designed 
by M.D. Burdick of this Bureau and heated by 
means of Globar elements The furnace is de- 
signed to accommodate six heating elements and 
= intended for use at temperatures up to 2,500° F 
Laodi> © 


used for calculating the values for the modulus of 


The quarter-point loading formulae 


rupture (ZR) and the modulus of elasticity 
VE) in bending were: 
» d6 yh 
Wik thd??? 
and 
WE Pash ha 


t Abd 
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where P, is the total applied load, L is the span, 
b is the breadth at the breaking point, d is the 
depth at the breaking pomt, A is the deflection 
for anv load 72, and a@ is one-quarter of the span L 

The interior of the furnace is shown in figure 2 
The specimen bar (4) was supported on free-titling 
refractory knife-edges (2), and the deflection at 
mid-span Was measured with a '\-in. diameter arti- 
ficial sapphire rod (C) and a calibrated gage read- 
ing to O.0OOL in. In order to eliminate any cor- 
rections due to the expansion or movement of the 
supporting knife-edges during the test, the gas 


and gage holder, outside of the furnace, were sup- 
ported with two sapphire rods (/?) resting on the 
specimen at the span pomts The four elements 
with which the furnace was heated for these tests 


were spaced symmetrically about the specimen, 


a 
E 





with two 


were obtamed 


and the temperatures 
chromel-alumet thermocoupl Ss; one located nbove 
and the other below the mid-span ol the specimen 


After having attamed approximately equilibrium 


conditions m= the furnace, the test bars were 
stressed in increments of 1,500 Tb ins every 2 min, 
and race deflection readings were taken every 


minute, When it was noted that the deflection 
following loading was practically constant by the 
end of the first minute, it was assumed that a 
permanent deformation, or “plastic flow’, had 


not occurred 
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IV. The System MgO-BeO-ZrO 
l. The System MgO-BeO, Fields of Stability 


he binary eutectic between MgO and BeO has 
been reported as beme between 1.670 ana 
1.855° C and 50 and 48 percent of BeO Ry Bas 
1s} \ more recent investigation [5] has shown 
the binarv euteetic to be at TS355° C and 48.2 


percent of 4) No compo ins or mixed crys 


tals were found 


2. The System BeO-ZrO., Fields of Stability 


In the system BeO-ZrO, a simple eutectic was 
first reported at about 2,180° C and 65 percent of 
ZrO, {14 The compound (BeO.2ZrO. and two 
eutecties at about 2,100° C and 65 and S85 percent 
of ZrO. were given in 1930 [15], but a vear later a 
single eutectic at about 75 percent ol ZrO. and 
2.240° © was reported [16 In the work at this 
Bureau [5] no evidence was found of a compound 
between BeO and ZrO 


3. The System MgO-ZrO 
(a) Fields of Stability 


hie phrase relations in this svstem are apparently 
complex, and many conflicting data have been 
ven in the literature Ruff and Ebert) [17 
reported two eutectics at about 26 and SO mole 
percent of ZrO) and 1,600° and 1,550° C and the 
NleO-ZrO. melting at about 2.150° C 


In acelition, thev indieated that a sertes of cubs 


COTTLPOULEN 


mixed ervstals exists between about 60 and 100 
mole percent of zirconia. ‘The first investigations 
of von Wartenburg and Werth produced similas 
results, but ina later re port [15 hes cave only 
one eutectic at about 50 mole percent ol zirconia 
and 2.160" Cooand no compound formation 
Ebert and Cohn IS corroborated these later 


data, and in 1939 Zhirnova [19] further defined the 
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simpli eutectic at 50 mole percent 79.3 werght 


percent) of ZrO, and 2,070° C 


In the present investigation no attempt was 


made to redetermine the binarv euteetic. but 


considerable number of X-ray studies of various 


mixtures were made in the Constitution and Micro 
structure Section of this Bureau. In general, the 
results of this work indicate that the conclusion 
ol Rut! and Ebert are correct in so far as the ¢ ube 
isomorphous series between 60 and LOO mole pet 
cent (S82 to 1L000,) of zirconia is concerned: but 
no compound of zirconia and magnesia wa 
identified whi li supports the work of Ebert aT 


Cohn 
(b) Ceramic Bodies 


Rieke and Ungewiss [20] made a number o 
ceramic bodies of MgO-ZrO, compositions, mainly 
for dielectric testing, and found that only on 
mixture (72.5) of MeO) in the range from 14 t 
100) percent of magnesia showed a porosity o 
more than O.L percent at maturity Shrinkag 
Values were rather high, from 20 to 28° percent 
The imperviousness of magnesia-zirconia comp 
sitions matured at high temperatures, 1,700 
1.S00° CC. was noted also in a previous study 
the thermal length changes of zirconia {9 

In this investigation only a few properties « 
selected body mixtures, ranging at about | 
pereent intervals from 3S to 92 percent of 7 
eotla were determined and these are riven 
table 2 All COT positlons between 50 and 9 
percent of ZrO. could be matured to Im pervied 
hess at temperatures between 1,700° and L.S00° ¢ 
The SS pereent of zirconia mixture was overtir 
at 1S00° C and was not completely matured 
1.750° © (0.150% absorption A complete stud 
was not carried out because of the relatively hig 
density of these bodies, although rather high valu 
29 000 to 30.000 Ib in were obtained for tl 
modulus of rupture in’ bending at  1,S800°) | 


982° C 
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4. The System MgO-BeO-ZrO 
(a) Fields of Stability 


The results of the fusion-point and softening 
range determinations and of the petrographic ATL 
\-rav examinations of some of the 61 mixtures 
nvestigated are summarized ino table 3 Only 
those considered essential to the establishment of 
the probable phiase relations are reported The 
the specimen that had been heated at or near th: 


when the liquidus was not reached. It must be 
emphasized that equilibrium may not have been 


attained and that the slow cooling (about 200° C 


per hour) in the furnace may have destroyed any 


evidence ol t hae phases present at the hiches 


temperatures 





a 
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recorded data for each mixture were obtained for 


estimated liquidus or at the highest temperature 


ystem MeO-ZrO 


Figcure > Is 2 graphic representation ol the SVs 
tem oas derived from the available data This 
system ts believed to be oa simple one with the 
ternary eutectic ator neal mole COPLPOstllon 


HNI:5SBi3Z and 1,669 i. _& The 


phase boundar 1s 


primary 
short-dash lines) are assumed 
to be essentially straight lines joining the various 
binary eutectics and the ternary eutecti No 
evidence ol compound formation was observed 
either with the fusion-point and softening-range 
determinations or with the petrographic and X-ray 
CNG CLOTS 

The intimately mixed oxides were used without 
presintering for the initial observations of thi 


melting characteristics in this svstem 
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Tari 3 Re nepoint ar oftening-range observation and of petrographic and ‘ mor ” 
compositions in the fe MgO0-BeO-ZrO, Continued 
{ 
Ma k raphie and X-1 ns (8 
mur \ i were 1 u 1, and the 4 
NI Weigt ! I i iv not b h brium at iN ‘ 
) BeO Zr N\lgeo Het) Zr 
( \ ( 
j si) aN _ s | i s he i rt X Nico, Bea 
il fu ‘ t Zr‘ 
j is Fu it “ | pal pl ur i iw i 
pole | ; yw il f i“ i ! 
' ' Fu I i SY Ma a, bo i, and i about se 
| ] ‘ i i il il 
near tl it 1 
s ’ Fu he " I ‘ 
it l 1! | 
i s Fu He l 
is s s s ij 2 > Vigo) 24Red) 
% } “ \ 
StS 12.8 4 I ‘ K blocky t s 
S . } s t} pt 
be ii “ i “ 
. t ! ‘ j 11 
? X ) Re Zro 
iM 5M \ Ve Reo ALD 
i Fu He 
i X 
oO. BeO ZrO 
s } i\lvo hea) 
tt \ { 
swell ( ? ) 
An attempt to identify the ternary eutectic mullite and boron carbide mortars to pass the 
solely by microscopic examination and subsequent No. $25 U.S. Standard Sieve. This finely divided 
nterpretation was only partially successful. The material was then pressed, shaped, and = tested 


most that could be done was to identify ah area 
the 
distinguishable 


n Which specimens examined contained no 


primary phases and which were 


composed mainly of a regular mosaic type struc- 


ture The area (fig. 3) so defined is bounded by 
the compositions 3M:5B:2Z, 1M:2B:1Z, 3M: 


5B:4Z, 1M :1B21Z, and 5M:5B:3Z. 
In order to determine the ternary eutectie more 
accurately, il number of compositions were pre- 


sintered at 1,.670° C for 30 min and ground in 
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the 


The results of tests on two of these presintered 


according to procedures given previously 


7) 


compositions are included at the end of table 
It was noted that the temperature at which fusion 


beg 


an for the presintered ternary eutectic mixture 


was 20° C lower than that for the unsintered 
muxture. 
An area of possible solid solution Is suggested 


for this system extending along the BeO-ZrQ, join 


from about QD percent to about SO pereent ol Be 
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and bounded by, but) exeluding, compositions The appearance of this exsolution texture is simila 
NIE:P4B IZ. IM :I1LOB:1Z \I:5B:1Z.3M:5B:4Z to that seen in specimens ol the \leO-Bet -ThO il 
and SALSSB:SZ. This area is shown on the dia system, but the texture is of much finer grain siz 


by a long-dash line In this area The appearance of lines for cubic ZrO, in the X-ra n 


ram enclosed 


there may be at least two phases containing limited patterns is not necessarily an indication of ternay 

nmounts of the end members mu solid solution eolid solution but rather the expected form ( th 

The ends of the line approximate the limits of ZrO, in the presence of magnesia |9, 19 oe 

solid sol iboal t\ Although thre ditferences inh the al 
vtive index values between MeO and BeQ (b) Ceramic Bodies ul 

oe AO _ — able, petrographic examina The only available literature information o ps 

tion failed to reveal any changes in these values ceramic bodies in this svstem /2] 22) Is concern 


mad only occasionally Were thre nalices of the cube 
form of ZrO sn) hou trv othae COMLPOsitlons 
INDS6BiSZ and \l:4B:4Z, however, ai ve 


mainly with modifying the thermal expansi: 


characteristics of zirconia by the add trom of sim: 


ry 
, ints ’ 1les 
ale fin) ti eNXsoiltittor tent Wis seer although Live ELOUT ol th thet ON af 
rvstals of lower refractive index than = zireonin Impervious bodies mav be made from com) 


vere too small to be sufficiently resolved by thy sitions in approximately one-half of the terna 


eroscope for complete optical identification diagram ares The compositions of these bod 
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Ficure 3 The system MeO-BeO-ZrO 
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lie to the right of the solid line shown in figure 3 


and could usually be matured between 1,450° and 
1.600° C, Table 4 gives the results of the matur- 
ing range studies and the values of porosity and 
shrinkage at maximum density. Those bodies 
that could be matured to a practically impervious 
condition over a range of 50° C or more have 


absorption values aver: 


ging about 0.02 percent 
and shrinkage values ranging from 15 to 21 
percent 


Faris 4 Vaturing range studies and compressive strength 


( | Pro 
Mok Weight Maturit 
sae les nye 
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Compressive strength values (table 4) ranged 
from 164,000 to 264,000 Ib in? for these porce- 
lains, and the outstanding body in these tests was 
IM:24B:1Z (264,000 Ib/in# The maturing 
range of this porcelain was considerable (1,450 
to 1,600° C), and its density (3.25 ¢ em*) was one of 
the lowest observed. The table vives some com- 
parative compression results for three body com- 
positions similarily fabricated but matured in a 
gas-fired furnace. 


esults of some compositions in the system MgO-Bet -ZrOvo 
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Nim Panty COMposi tions were sS¢ lected to repre- 


sent the relative resistance to thermal shock. and 


strengths im bending at room and at elevated 


temperature of the refactor porcelains in the 
svstem Table 5 gives the data for these proper- 


ties and also includes some values of the modulus 


of rupture and Young's modulus for a few bodies 
tested at 1,900° F(1,038° C) and 2,000° F(1,093° C 
As expected, the high-bervilia bodies were the 
most resistant to thermal quenching, and the 
high-zirconia bodies either did not survive the 
heat-shock test or were very easily broken by 
linger pressure alone after testing. Generally 


the rupture values of the high-bervllia bodies were 


shehtly higher after sho« king than those obtained 
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for unshocked specimens. The values for the 
modulus of rupture in bending at room tempera- 
ture for the high-bervllia bodies ranged from 
13,000 to 20,000 Ib in, and from 23,000 to 39,000 
Ib in for the high-zirconia bodies Modulus of 5 
rupture and modulus of elasticity (Young’s modu - 
lus) values in bending at 1.800° F for the high 
bervllia bodies ranged from 13,000 to 18.000 
lh in. and 38,000,000 to 34,000,000 Ib in’, respe 
tively, and for the high-zireonia bodies the values ) 
were from 21,000) to 35,000) Ibin’ and from v 
31,000,000 to 15,000,000 Ib/in2, respectively. It 
is Interesting to note that, at 1.800° F, the Youne’s 
modult values decrease as the moduli of rupture 
values Increase 
ri 
gtl npositions in the system MgO-BeO-ZrO \ 
| 
y g’sM ol 
tt ober! ; ' [ 
\. 1} - peratur 
PLIALL | 
2) oo 7 l 
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Kn . | 
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wT . min | 
' " nan 3 | 180 
; ‘ 8.4 1 
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V. The System MgO BeO ThO 
1. The System MgO ThO., Fields of Stability 


Von Wartenburg and Prophet [11] deseribe the 
system MgQO-ThO, as having no melting below 


> 500° C 


2. The System BeO ThO., Fields of Stability 


\ simple system of one eutectic at about 20 


percent of BeO and 2.200° C 


Was reported hy von 


3. The System MgO BeO ThO 
(a) Fields of Stability 


Wartenbure and coworkers 


The results of the fusion-point and softening- 


range determinations and of the petrographic and 


X-ray examinations of some of the mixtures 
investigated are summarized in table 6. Only 
those considered essential to the establishment 
of the probable phase relations are reported 


The recorded data for each mixture were obtained 


for the specimen that had been heated at or near 


BLE ig mi-point and oftening-range obse 
”" asil , nm the 
( 
Ma 
mul Effect 
NI \ t" 
All te 1 
te 
ure 
let) hk ) ‘) Nieo) Reo Tho 
f 
’ wsD SZ k t ina 
| Lat 1.832 
y iz Hoy S71 I n bega 
need fu 
2 12 Hs s Fu nt ma 
u it 1,87 
s “44 23. 22 ; Fusion t ina “ 
1 al 
2s H 4 l 2 “MW Fu it in at al 
1S, mo furtt 
up to 1,900 
‘ iu wT) sid Fu 1 began at 
| i si 
} ! ‘ sy Fusion t wn a s 
pleted } 
i ( SIS I it i 
See footnotes at end of table. 
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the estimated liquidus or at the highest temper 


ature when the liquidus was not reached. It will 
bear repeating that equilibrium may not have 
been attained and that the slow cooling (about 
200° C per hour) in the furnace may have de- 


stroved any evidence of the phases present at 
the higher temperatures 

Figure 4 is a graphic representation of the svstem 
as derived from the available data. This system 
is believed to be a simple one with the ternary 


eutectic at or near mole composition SME LOBITT 


and 1,797 5° °C. The primary phase bound- 
aries (short dash lines) are assumed to be essen- 
tially straight lines joining the various binary 


eutectics and the ternary eutectic. The boundary 
line between the ternary eutectic and the MgQ- 
ThO 


mating the results of this investigation. 


eutectic is shown partially drawn, approxi- 
No evVvi- 
dence of compound formation was observed either 
with the fusion-point and softening-range deter- 
minations or with petrographic and N-ray exam- 


Ihatlons 


or and of petrographic and j ne 
fem MeO-BeO-ThO. 
Resu | 1 X-ray exat 
“I er ‘ ‘ not quenched and the pl t 
( t ed l t in the x 
i ! yu 1 
B the | ry pl 
ul | wl ! = be 
Irs «as ¢ ition from bery xX weak MgO 
BeO 1 rho 
s lar to ce 2MgO:7BeO eT ThOs, except that Ooh 
the X-ray shows medium stro rho 
s rt , 2MyvO-7BeO-tThOs, except for 65s 
! beryilia by px rapehiie 
Kye i sith i i 
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Beryllia ils conta both i rand th i 
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sth i t rt ear the ternary eute X 
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Ma " he prima phase lhorta i i 
‘ v la Lo t 
il 
ippea be the | | phase I 
" xsolut beryl si } 
ipl j } ym { 
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All compositions for the fusion-point and 


softening-range determinations reported in this 


system were presintered at 1,750° © for 30> min 
in mullite and boron carbide mortars 


This 


shaped, 


and cround 
to Piss thre No 


finely divided material was then pressed 


25 (| S. Standard Sieve 
ana te sted necording to thre procedures previously 
dese! bed 

ternary 


The location and temperature of the 


eutectic of this svstem were based almost entirely 
upon the data obtamed from the fusion pom and 
soltenmg-range observations, because petrographic 


examinations failed to reveal even an area sul 
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if 


rounding the probable ternary 


found for the zirconia system 


of the specimen cOMpostllons in all sections of tl 


diagram 


tvpe structure as seen 


A ana 


prac tic 


=] owed both 


B oof figure 5 


In addition, mar 


2 mosaic and i 


NeQ-BeOd-ThoO Continued 
K " ! iX ‘ n 
( t | be t i ul \ 
the 
i 
ait? \ 
pt } , 7) ‘ 
‘ x 
i ippears to be the | ry phase t ‘ 
I ie t he ul at 
“w i " i? 
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ho 4 k MeO 
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eutectie as W 


in the photomicrograp! 


lly always occurred 


dendritic 


nto oan 


struct 


Wnheasure 





depth Ith those surfaces exposed lo th furna 
atmosphere, and the mosaic structure usua 
constituted the remainder of the thin-secti 


specimen The photomict wraphs also rive 
indication of the relative difficulty in distinguis 
ing small or large crystals of either magnesia 


bervilia without polarized light 
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Figure 4 The system MgO0-BeO-ThoO, 


/ 


figure gives the locatio fthe primary phase boundaries hort-dast uL tu i 
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the system MeO-BeO-ThoO 
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Because of the apparent highly viscous nature 
of compositions in this system at high tempera- 
tures the re still retnaritis considerable doubt 
whether the ternary eutectie may be located more 
accurately than in the area (fig. 4) bounded by the 


TVML:LIOB:1T. SM:i1OB:1T. | and 


COMLPOstlLons 


ONL:9B:27T The  fusion-point and = softening- 
range tests however, indicate it to be near mole 
cCOMposi tion SM:1IOBIUT and 1.797 S* & 


An area of possibl solid solution is suggested 
for this system to be bounded by and to include 
COMpPoOsttLons 12\I:78B:1T 6M:42B:1T, 2M: 
lOB:LT, and 3M:28B:47T This area is shown on 
thre diagram enclosed by a long-dash line The 
X-ray examinations failed to reveal any shift in the 
l-spacing values, and the petrographic examina 
tions failed to reveal uns differences in refractive 
index values ol the component oxides, although 
a definite exsolution texture (photomicrograph B 
all compositions within 


of fig. 5) was noted in 
the area Examination of the softening-range 
observations (not all of which are included in 
table 6) shows that the liquidus temperatures of 
the compositions located in the central portion 
of the diagram are between 1,814° and 1,840° C 
ana that, ut the proposed solid solution boundary, 
there occurs a large increase in these temperatures. 
Within the solid solution area, the mixtures of 
the compositions 12M:78SB:1T, 6\L:42B:1T, 
2M: 15B:1T, 2M:ilOB: LT, and 3M:28B:4T show 
only very limited melting in the temperature 
range from 1,870° to 1,915° C, whereas those 
mixtures just outside of the area, 30M :120B:1T, 
I2\E:42B:1T, 6M:30B:1T, 6M:20B:1T, and 
2N1:7B:1T, are completely melted in the tem- 
perature range from 1,812° to 1,842° C 

It is believed that the very viscous nature at 
high temperatures of the mixtures in this region 
has obscured to an appreciable extent the expected 
increase in the solidus temperatures associated 
with the solid solution phenomena. All of the 
available data seems to indicate that the suggested 
solid solution area is one that exists only at high 
temperatures, and that the solid solution is not 
present after slow cooling, which apparently 
destroys all evidence of the existence in the speci- 


mens of the high-temperature form 
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(b) Ceramic Bodies 


It is possible to produce impervious porcelains 
from compositions in approximately one-half of 
the beryllia field of this system These bodies, 
whose compositions are confined to the solid 
solution area, could usually be matured between 
1,500° and 1,750° C Table 7 gives the results of 
the maturing range studies and the values of 
absorption and shrinkage at maximum density 
Those porcelains that could be matured to a 
practically impervious condition over a range of 
50° C or more have shrinkage values ranging 
from IS to 21 percent and absorption values 
averaging about 0.02 percent 

Compressive strength values (table 7) ranged 
from 165,000 to 265,000 Ib in’ for these bodies 
The outstanding porcelain from the standpoint 
of compressive strength was 3Mi28B:4T (265.700 
Ib in.*), but its maturing range is relatively short 
LOlecom 


1.600 to 1.650° (C) and its density 


was the highest observed 

It had been found previously [5] that thoria 
containing bodies were extremely brittle and 
usually possessed rather poor heat-shock resist 
ance. For practical purposes it was, therefore 
considered sufficient) to determine the relative 
thermal-shock resistance and the = strength iy 
bending at room and at elevated temperatures fo 
only a few bocly COMpositions selected to represent 
the maturable area of the system. 

As had been expected, none of the bodies test: 
(table 8) were satisfactory in their resistance ¢ 
thermal shocking (all were either cracked 
could easily be broken after testing) nor coul 
their strengths in bending at room temperatut 
9000 to 14,000 lb/in2) be considered as mor 
than fair to poor. Modulus of rupture ar 
modulus of elasticity (Young's modulus) valu 
in bending at 1,800° F (982° C 
13,000 to 16,000 Ib/in’ and from 40,000,000 
$2,000,000 Ib in The higher valu 
for the tests at 1,800° F as compared to thos 


ranged fro 
“. respectively . 


obtained for the room-temperature tests ma 
have been due to the incomplete maturing of : 
but one of the test bars. The values obtain 
for body composition 2M:42B:1T, which w 
completely matured, were slightly lower for t! 
1.S800° F test. 
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VI. Conclusions 


No outstanding porcelain-ty pe bodies were 
found in either of the systems studied, which com- 
bine all of the desirable high-temperature proper- 


ties required for their use as components of heat 
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engines, although a number of bodies are now 


high- 
temperature services Where thermal-shock condi- 


available for various high-strength and 
tions are not of prime importance 
The limited of the phase equilibrium 


relations suggests that the svstem MgO-BeO-ZrO. 


study 
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be considered as a simple one with the ternary 
at composition 4\MgO:5BeO 3ZrO, and 


1.669 a S| 


eutect& 
No evidence of binary or ternary 
compound formation was observed and an area 
of solid solution is suggested to extend along the 
BeO-ZrO, join 
matured at temperatures ranging from 


Ceramic bodies in this system 
may be 
1.450° to 1,600° C in approximately one-half of 
the ternary diagram area The values of the 
physical properties of these bodies are summarized 
in table 9 

The summary table of these physical properties 
is included in order that the veneral properties ol 


these porcelain bodies may be readily compared 


with the properties of other refractory bodies, 
such as are given in table 1 
Pantie Summary table of tt yeneral physical prope 
if ; } j po hod 
\led)-Z7r¢ f VieQ)-BReOQ-Zro VieO-BeO-Tho 
' o ~“ ; J 
i s iu 
\l i , e 
wre f ‘ 
( 
mi } 
i 
i 
Some limited data were determined for a few 


porcelains in the binary system MgO-ZrO,, and 


these are included in the summary table 
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The phase relations and the location of the 
primary phase boundaries in the system MgQ- 
BeO-ThO, are similar to those in the previous 
system, although this is a distinetly high-tem 
perature one with the ternary eutectic at about 
5MgO:10BeO:1ThO, and 
1797 a € No evidence of compound for- 


mole — composition 
mation was observed; and an area of solid solution 
is suggested, bounded by an S-shaped curve ex- 
tending from about 25 percent of MgO on the 
MgO-BeO join along the line of 20 percent of 
MgO to about 70 percent of ThO, on the BeO- 
ThO, join. Ceramic bodies in this system may 
be matured at temperatures ranging from 1,500 
to 1,750 C 


bervilia field of the system bounded by the samy 


in approximately one-half of the 
line describing the suggested area of solid solution 
The range of values determined for the various 
physical properties of these matured bodies 
summarized in table 9. 

It is of interest to note that the maturing ranges 
of the porcelains of each of the systems are such 
that all of those studied in the system MgO 
BeO-ZrO,, some of those studied in the system 
MgO0-BeO-ThO,, and none of those studied 
MgQO-ZrO. can be 


kilns currently 


the system matured im thi 


commercial being used in) the 


porcelain industry 


MeMurdic 


Sureau’s Microstructure 


The authors are grateful to H. M 
and Ke 


Section for the preparation and interpretation of 


Golovato of the 


the X-ray powder patterns and to P. J. Yavorsky 
formerly of the Bureau staff, for a portion of thi 
preliminary investigation of the system MgO 
BeO-ZrO,, 
VII. References 
1} Gas Turbines and Jet Propulsion, NBS Letter ¢ 


cular LOCS72 (Auge. 6, 1947), and semiannua 
2} R. F. Geller and M. D. Burdick, NACA ARR6D2 


3} J. R. Bressman, J. Am. Ceram. Soe. 30, 145 (1047 
1} Ht. C. Cross and W. F. Simmer OSRD No. 4717 
ser. No. M-477, NDRC, OSDR, War Meta 


R. F. Geller, P. J. Yavorskyv, B. L. Steirmat 
\. S. Creamer, J. Resear 
RPI703 
6) F. P. Hall and H. Insley, J. Am. Cerar Soc, 30, 
11) Part IL (1947 


Journal of Research 














dards Miscellaneou 
10, 1947 


Yavorsky, J 


National Bureau of 
tion MIIS3 


oO) RI Cseller ar 


“tear 
a 


iP.d 


S7 (1943) RP1662 
) \pparat leseribed ty 2.'4 Harr J. Am. Cera 
Soc. 8, 774 (19025 IS 
Il ‘ Wa t and | Pr ) Z. a ! ( 
208, 360 (1932 y 
©. Ruff, I bobo ind Krawe Za i 
( 213, 333 (1933 20) 
s ol \\ i J Re and | Zara 
>] 
hs. & vy. a ( 230, 257 (1037 
1] O. Ruff, | bt ind k.. Stepha 4. a i 
( 185, 22 at) 


Refractory Oxide Porcelains 





von Wartenburg and 


Chem. 196, 374 (1931 
180, 


©. Ruff and F. Ebert, Z. anorg. allaem. Chet 
19 (1929); 87, LOS (1914): 88, 39S (1914 ind 133, 
193 (1024 

KF. Kbert and k. Cohn, Z. a i Che 713, 
j21 (19838 

N. A. Zhirnova, J. Applied Che SSR) 12, 127 
1a30 

R. Rieke and A. I 1h By ra ( 
17, 237 (1936 

| Rvs ( | s 4 t No. 1.060.000 (A 
7 Los4 a i Cer i Pate No 925.620 
Apr. 8S, 1930 


Cc. E. Curtis, J. Am. Cera 


»> 


WasHiIncron, April 1949 


447 


























Fle 





U. S. Department of Commerce Research Paper RP2035 
National Bureau of Standards Volume 43, November 1949 


Part of the Journal of Research of the National Bureau of Standards 





Laboratory Flow Tests of Fixed Spray Nozzles with 
Hydrocarbons and with Air 


By M. R. Shafer and H. L. Bovey 


Phe metering characteristics of fixed spray nozzles of the t vy ised in some turbo-jet 
nvines have been investigated Some of the nozzles supplied by the Navy Department 


ntained burrs, metal particles, and improper machining, which caused erratic fluid meter- 





ny \fter being reconditioned, a group of 26 nozzles was flow-tested with five different 
fluids to determine the effects of fluid density, viscosity, and supply pressure upon the rate 
of discharge of the nozzles The results indicate that it is impracticable to correct for 
differences in the physical properties of the test fluid \ comparative method of flow- 
testing fixed nozzles with air is described Although this method leaves much to be desired, 
it appears useful for safe and rapid sizing of nozzles to within ~ 3 pereent of their actual flow 
I. Introduction months In order to study the effects of the 


" : densitv and viscositv of the test fluid, they have 
For several vears the Bureau of Aeronautics, : 
1) f \ l l been tested repeatedly over the pressure range ol 
epartment of the avv, has sponsored at the " ‘ 4 » : 
| - | , » to 250 Ib in with the following liquuid hvadro 
National Bureau of Standards a program of testing ; Pp 

' carbons: Varsol, pure n-heptane, Apeo 467 oil, a 
and research on devices for handling and metering 
ls f ft \ t pl ft commercial mixture of isooctanes, and Soltrol 100 
uels for airera A recent phase of this program 
Is concerned with the flow characteristies of fixed In addition, an attempt has been mace to 
spray nozzles of the type that have been used in develop a method using air, instead of a flammable 
SOT turbo-jet engines. Depending upon the sup- hydrocarbon, as the test fluid, Phis Fepors pre- 
plier, these may be designated as Monarch nozzles, sents the results that have been obtained to date 


Hago nozzles, ete., and for the present purpose all with the six different test media 
are essentially alike 

For the present work, three sets of sixty nozzles II. Description of Nozzles 
each were procured. Each set bore a color desig- 


hation indicating that all of the nozzles of each As shown mn igure I, the nozzles CONSISL essen- 





set hed been matched in few to within 2.5 per- tially of a body, an insert, and a strainer. hes 
cent at a pressure of 100 Ib in are designated by Navy Parts List No. 14G520-4. 

As a first step, all nozzles were flow-tested as After passing through the strainer, the fuel is 
received with Varsol at five different’ pressures directed by tangential slits in the insert into the 

It was found that a large proportion exhibited swirl chamber formed between the end of the 
significant changes in flow characteristics before insert and the body. The kinetie energy of the 
and after being tested at a pressure of 250 Ibiin- fuel in the swirl chamber is effective in atomizing 
Subsequent disassembly showed that these changes it as it escapes through the orifice in the body 
were caused by burrs, metal particles, and poor The dimensions and relative locations of the tan- 
Inachining gential slits, the swirl chamber, and the orifice 
From the lot of 180 nozzles, 26 were selected, determine the pressure-flow characteristics, the 
cleaned, reconditioned, and assembled for use in spray angle and distribution, and the drizzle point 

the remainder of the tests. The flow characteris- of the nozzle. The latter may be defined as the 

ties of the 26 have remained constant over several lowest fuel pressure at which the nozzle produces 
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av, and below which it discharges laren: drops 


ora stream of fuel 


III. Flow Tests with Varsol of 180 Nozzles 
as Received 


Figure 2 isa diagram of thi upparatus used for 
flow testing the three sets of 60 nozzles each. as 
they were received. The Varsol was circulated 
bya pump through a low-pressure cirenuit from a 
storage tank, thro igh a heat Ne hanes and back 
to the tank ata constant rate of about 200 eal bin 
Fuel to the nozzle was bled from this line, sunne 
passed through a Rotameter to a second pump 
having UpPproapriate Valves ina by Puss ana ith th 
discharge line From this pump the Varsol Prussed 
through a LO-micron filter to thy fitting bearing 
the nozzle Fuel pressure Was measured at) this 
hitting 

The Rotameter was calibrated with the fluid 
ised for the tests and a thr temperature of the 
tests. Observed values of flow ar believed ae 

irate to within at least 0.25 percent The 
pressure gages were also calibrated at) intervals 


titiad gre beleved to bye necurate to within () .) 


Ih oan over thre range trom YS to 250 |b in 
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The 180 nozzles were tested as received wit! 


Varsol having a kinematic viscosity of | 102 cent) 
stokes and a specific gravity of O.779 at the test 
temperature of SO° F, The flow. of each wa 
measured at the following pressures and in th 
order stated: 98, 150, 250, 50, 98. and 5 Tb in 


race The results 


except for the initial valu 
at OS [bin are shown in figures 3. 4 and 5 fo 
the sets of 60 nozzles color coded red Creech, ane 
purple, respectively The flows shown at Os 
Ih in were observed after the nozzle had bee 


~ tbjected once to thr higher pressures 
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Specifications for these nozzles state that the 
flow at 250 Ib in’ shall be in the range from 154 
to 160 percent of the flow at 100 Tb in’. Correct- 
ing these limits to the pressure of 9S Ib in’ used 
in the present tests, the limits become 156 and 
162 percent of the flow at the test pressure. Speci- 
fied flow limits for these nozzles at 9S and 250 
lb in? are indicated by dashed lines 

It will be noted that some of the nozzles failed 
to flow within the specified limits, particularly the 


» 


red group (fig. 3) at 250 Ib in It is also ap- 
parent that nozzles that are matehed in flow at 
one pressure frequently are unmatehed at othe: 
pressures, 

Many of the nozzles flowed differently at YS 
Ib in. before and after being subjected to higher 
pressures As will be seen in figure 6, this change 
continued for many pressure eveles with certain 
nozzles. The zero lines represent the flows at 
YS |b in shown in the previous three figures 
Flows at this pressure changed by more than | 
pereent for 46 of the nozzles as a result of subject 
ing them to pressure eveling. This indieated 
that something must have been moved about 


Within the nozzles during the tests 
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IV. Reconditioning of 26 Nozzles 


) 


Following the above tests 26 of the nozzles 


were taken apart snd the following defeets were 


found in all to varving 


extents 


There were meta! chips resembling filings in 
the tangential slits and swirl chambers; 2. Many 
, 


of the seals between inserts and bodies seemed 


\lost of the 


finished on the sealir end, and the edges of tha 


mperfect nserts were not properly 


slits were rageed t. Large burrs left in cutting 
tha lits remained attached and caused partial 
blocking of the swirl chambers 


In the enlarged photograph shown in figure 7 


such burrs are visible ms ure thre rough surtace o 
thr Insert ana thy ragved edeves of thr slits 
Tha burrs nthe chip Vel removed hie each 


th its original insert 
No. sienifieant 


occurred subseq rently “> that t seems salt to 


nozzle was reassembled w 
nha strainer change in flow has 
state that nitial changes were due to movements 
of burrs ind chips by the test fliricl lt is obvious 


Presence ot si hy fore rn particles enbnotl 


V. Flow Tests of 26 Reconditioned Nozzles 
with Hydrocarbons 


Phe 26 reconditioned nozzles. renumbered in the 


orc ol mereasig 


flow enpaeirty with Varsol ut OS 
Ih in.-, Were next flow-tested at pressures of 9S and 
250 tb ine with five a hydrocarbons having 


the follow propert 


Test fluid No. 1, designated as Varsol. is a ven 


erat utiiityv sorvent stocked as a storeroom tem 


at this Bureau It is similar to cleaner’s naphtha 











hiatre 7 D f no 
Test thiids No. 2, 4, and 5 were supplied through | 
the courtesy of the Phillips Petroleum Co il ] 
test fluid No through the courtesy of tl 
Anderson-Prichard Oil Corp ) 
Fluids No. 1.2. and 3 were selected because ea 
s thought to be under consideration as a standa 


I] lltis AY 


band 5 were selected because thev have the san 


fliid for te sting yet engine nuNXiliaries 


densitv. but widely different viscosities Of 1] 
ve heptane nod the mixture of tsooetanes h z 
properties approxtmatin those of aviation gas 


line, Whereas Varsol and Apeo resemble keroses 
Thiaore closely 

The results obtamed with 26 nozzles at pressul 
of OS and 250 Ib in for each of the five Hau 
hydrocarbons are presented in table | 

igure Sis prese nted to show the re prod iibil 
of the flow measurements and of the nozzles s 
sequent to reconditioning Kach point represe 
the deviation of one observation made with Va 
from the average of all observations made w 


the same nozzle, fluid, and test) pressure | 


mnaXximum deviation does not exceed 0.5 per 
from 


Thus : 


as well as the Iheasiurenvet 


and only 5 of 100 observations deviate 
mean by more than 0.25 percent 
nozzles themselves 
uppear satisfactory for present purposes 
effected by 


The changes reconditioning 


shown in figure 0 As will be se 


nozzle s ure 


Journal of Resear: Fl 














Pane | Flow « q cart 
~ ~ } . 

' 
> —_—_ — —————————————— a- —--— + - $$ $s —_ _ — 
" | . 

: j : 
f . . ‘ . °° $° 
: L ~- —_ . ; dus » S21 —S 4 é ° 
| . - 
| - 7 4 ; 3 o*?? o 9 
5 5} % eC 6 a etee8 G 
z 
=e = es? 
m 2 3 - a ~ 24 . . 
NOZZ NUMBER a 
| RE S \ j } | iS 
] te \ ‘ 
) a } A 7) a — = -" - 
. MBER 
a 7 | 10 Flo a , nana g 
4 @. Varsol; 5 
Nee the changes in flow ranged from —-S to —S percent, 
* and only three of the 26 nozzles changed less than 
| percent 
: . The data of table 1 for Varsol, Apeco, nnd hie » 
= 6 tane at a pressure of 9S Tb in’ are compared im 
| | | figure 10. The nozzles were numbered arbitrarily 
i A On elles ill in the order of increasing flow with Varsol, so that 
NOZZLE NUMBER 
avReEe 9. ¢ : are | , ning a relatively smooth curve is obtained for this thuiid 


Varsol at 9S psig lLlowever the corresponding curves for the other 
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two fluids are not smooth, indicating that nozzles 
matched for one fluid are not necessarily matched 
for fluids having different properties, and showing 
the need for caution in selecting a standard test 
thiviel 

To further emphasize this pom, consider the 
With 
an average flow of 41.57 


\\ ith 


average flow 


example furnished by nozzles 6, 7,8, and 9 
Varsol these 
) 


iD hie nnd i total spread ol O25 percent 


nozzle show 


heptane the order is different, the 
| pereent lower, anc the spread is 3.7 


With peo the 


the average flow is 1.4 


pereent 
flows are in a still different ordet 
percent lower, and the 
There are 
similar examples in the data of figure 10 

From the data of table 1, the 


0 Tb mn 


spread is 3.7) pereent many other 


ratio of the flow 
nt can be calculated 
For fluids 


this observed ratio is as follows 


to that at OS Ib in 
for each nozzle and for each test thud 
No. 1, 2. and 

Varsol 


j 


aving LO0.2 percent; 


ranging from 157 to 164 percent, aver- 


ranging from 156 


it prtaarne 


to 161 percent, averaging LOS.7 percent 


>and (peo 


from L600 to 167 percent, averaging 163.2 


\s a further wbustration., the ratios are plotted 


in fier In the form of a frequeney curve with 
, 
3 
we * 
7 
+ vw — 
* . © 
‘ n 
J 
a“ ~ ~ ; 
200 
hy ' j j 
Os 
= 


each point showing the number of nozzles having 


from. these 
fli 
ratio can be specified 


fluids No 


chosen to show the 


n given ratio It will be apparent 


results that the characteristics of the test 
must be known before the 
(As stated previously fand 5 have 


the same density and wer 
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effects of viscosity on nozzle performanc: 
12 shows how the results obtained with Soltrol 


varied from those obtained with the mixture of 





imooctanes It is at once apparent that the effect 
of viscosity varies in both magnitude and direction 
Henee it 


means of correcting for 


from nozzle to nozzle does not seen ) 
possible to develop any 


the VISCOSILV Of the test fluid 
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VI. Flow Tests With Air 


Karly in this nozzle test program wuppeare 


desirabl to attempt the development ol equ 


ment sine al instead ola flammable hvadroe arbeo 


us oa test thiid Such development seemed 


logical extension ol thre previous successful evol 


tion of the Navy Orifice Comparator, in which a 


s used lo flow test metering jets of atreratt en 


byugee tors with erentel speed necurney nnd safe 


than could be attained by other methods 


In applying the method to spray nozzles, 1! 
apparatus shown diagrammatically ino figure 


Was Investigated sriefly, air) compressed 


90 Iban. or more is passed through a pressu 


regulator and a filter to a test fixture consistin 
essentially of two chambers, each about 1'y in 

diameter by 4 in. in length, separated by a sma 
orifice or bleed The second chamber serves as 


mounting for the nozzle, through which all of tl 


regulated air es« upes to the atmosphere Pressu 
taps and manometers provide for observing tl 
pressure in the second chamber and the drop 
pressure between the two chambers 

In operation thre presstire in the second chamb 
is held constant, which means that the drop 
pressure across the nozzle Is also constant, and 
drop In pressure across the bleed between the t\ 
chambers ts observed. The latter is determined | 
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Figure 13. Schematic diagram of nozzle air lest, 


the volume rate of flow between chambers, which, 


in turn, is a function of the flow characteristics of 


the test nozzle 

Obviously the method is comparative rather 
than absolute, so that a set of nozzles calibrated by 
some other method is required for the calibration 
of the equipment using air. For nozzles having a 
flow capacity in the range 40 to 50 lb/hr at 100 
Ib in, it has been found by experiment that the 
air-test method gives best results when the bleed 
between the two chambers consists of a single hole 
made by a No. 76 drill (diameter=0.019 in.), and 
when the pressure in the second chamber is from 
1.2 to 2 times the pressure drop across the bleed, 

In developing the air-test method, the 26 
nozzles mentioned previously were tested with 
various constant pressures in the second chamber 
throughout the range from 10 in. of water to 50 
in. of mercury. The best) correlation between 
results with air and with liquid hydrocarbons is 
obtained in the range 25 to 35 in. of mercury. 
An example of this correlation is given in figure 
14. in which the pressure in the second chamber 
P,) was 30 in. of mercury, and the observed drop 
in pressure across the bleed (2?,—P,) is plotted 
against the observed flow of Varsol at a pressure 
of 98 Ib in’. 


observed points seems to be a straight line, from 


The best smooth curve through the 


which the maximum deviation is less than 2. per- 
cent and the average deviation ts less than | 
percent. The sensitivity of this air apparatus 
was about | in. of mercury per pound of fuel dis- 
charged through the nozzle per hour. 

The dashed lines in figure 14 are the limits of 
flow specified for nozzles coded red, green, and 
purple. On the basis of the results shown in this 
figure it might be concluded that the air-test 
method was satisfactory for fiow-testing spray 
nozzles. If this were true, it is certainly to be 
preferred from the standpoint of speed and safety. 

As already stated, figure 14 compares the 
results obtained with air and those obtained with 
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Varsol. 


of results with air and Apeo, and figure 16 shows 


Similarly figure 15 shows a comparison 


a comparison of those with air and #-heptane, all 
data for liquid fuels being for a pressure of 9S 
lb in’. Deviations from the straight: lines are 
all within +3 percent except for one nozzle in 
figure 16. Considering that the nozzles do not 
perform consistently with the three liquid test 
fluids, as was shown in figure 10, the correlation 
of the results with air and those with liquids is 


better than might have been expected. 
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In its present state, the air-test method would 
seem to have considerable merit for preliminary 
particularly in’ production and 


tests of nozzles, 


after overhaul. It would certainly be valuable in 
reyecting nozzles that flow so far from the design 
value that tests with liquid would be a waste of 
time, and in grouping nozzles within reasonably 
also be useful in 


limits. It might 


the matching of nozzle bodies and inserts to get 


narrow flow 
desired performance 

As a matter of fact, the uncertainties in the 
performance of nozzles determined with air are 
probably little, if any, greater than the present 
uncertainties in the actual performance of nozzles 
in engines. Consequently the further develop- 
ment of the air-test method will be carried along as 
rapidly as the resolution of the over-all problem of 


nozzle performance secs to warrant 


VII. Discussion and Conclusions 


The nozzles used in these tests were received 
late in 1947, and are believed typical of nozzles 
that The 


presence of chips, burrs, and inadequate finish on 


of this type in production at time 
certain parts were pointed out months ago, and 
it seems probable that such easily remedied faults 
mav no longer be a matter of concern. 

A fixed 
tions, namely the production of a spray and the 


nozzle performs two important fune- 


metering of fuel. This report deals only with fuel 


metering In nozzles of the type tested, this is 
accomplished primarily at two locations within the 
These ure oat the tangential slits 


nozzle foul 


that 
orifice in the nozzle body, which is in series with 


operate in parallel, and at the discharge 


the four slits It is not therefore, 


that the simple law of discharge through an orifice 


surprising, 


does not hold exactly for the nozzle as a whole 
the relation between pressure and 
flow for nozzle iter 
that of another, and the effect of a given change in 


(As examples, 
one rans somewhat from 
Viscosity of the test flriicl Than bye Ih one direction 
for one 


another 


Ixamined from this point of view, it is surprising 


that the results obtained with six different. test 


fluids including air agree as well as they do 


Perhaps the most obvious conclusion from these 


studies involves the choice of a fluid suitable for 


testing fuel nozzles and other engine auniliaries 
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nozzle and in the opposite direction for 


It is believed that a fluid chosen as a standard for 
this purpose should have as many of the following 
characteristics as may be found attainable: 

(a) It should be cheap and free from hight, 
strategic ingredients; (b) It should be availabk 
in quantity from many sources, and its physica! 
properties should be readily reproducible from a 

safe, 
effects 


as engine fuels on packings, diaphragms, ete 


production standpoint; (¢) It should be 


noncorrosive, and should have the same 
(d) Its physical properties should be the same as 
those of engine fuels and should not change with 


light Why Nn 


exceptions are made to the above requirements, it 


use by evaporation of ends; (©) 
must be possible to interpret data obtained with 
the test fluid in terms of performance in engines 

If the test fluid differs in density and or vis 
cosity. from the fuel used in’ the engine, the 
present results show that: 

(a) The actual rate of fuel delivery to the engin 
at any pressure may differ by several percent from 


Nozzles 


matehed for the test fluid at a particular pressu 


the rate predicted from the test data; (b 


may not be matched at this or any other pressur 
in the engine; and (c¢) it will not be practicable to 
develop corrections which are generally applicabl: 
for differences between physical properties of thy 
test fluid and the engine fuel. 

Both density and VISCOsILV of the fuel being 
metered by a nozzle are important in determining 
Bot! 
those properties change considerably with ten 
Thus it that 
might be by giving mor! 


its rate of discharge ata viven pressure. 


perature, seems likely turbo-j: 


operation improved 
attention to the temperature of the fuel enterins 
the individual nozzles, and to possible method 
for its control. 

Much additional thought and experimentatio 
are warranted in the development of nozzle tes 
data that will be truly indicative of the perforn 
anee of nozzles in operating engines Kv 
though some fluid is selected as a standard for t 
remains the development 


purposes, there 


satisfactory test methods and the evolution o 


significant test specifications. An obvious firs 


step toward the latter is the determination 
tolerances in flow that can be allowed in engin 
and more particularly of the unavoidable dith 
enees encountered in operating engines 

Thus the over-all problem of determining, | 


means of bench tests, whether a given set of spr 
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nozzles will give optimum performance ino an 


engine is highly complex. 


will be no question that they will perform as well 
as the original set 

sO long as spray nozzles are used in engines, 
that 
mentioned problems will pay dividends in’ im- 


there is no doubt solution of the afore- 
proved engine performance, particularly at high 
altitudes 


that 


Employment of spray nozzles requires 
level of fuel 


which in turn involves mechanical difficulties with 


the general pressure be high, 
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Equally difficult is the 
testing of a set of replacement nozzles so there 


fuel lines, manifolds, and seals, and in- 


creases the fire hazard in case a fuel line is broken 


pumps, 


It. therefore, seems legitimate to raise a question 
as to the relative merits of expending research 
effort on the further development of the various 
components of high-pressure fuel systems, or on 
the development of combustion chambers that 
will function with low-pressure fuel and without 


spray nozzles 


Wasnincron, January 10, 1949, 
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Mechanical Production of Very Thin Oscillator Plates 
By Leland T. Sogn and Walter J. Howard 


To produce quartz oscillator plates thinner than 0.005 inch, it has been necessary to 
modify the usual lapping methods and machinery. The thickness at which the erystal 
carrier (nest) stretches or buckles is the limiting factor in the conventional apparatus. The 
modifications described in this report consist of replacements for the conventional top 
lapping plate with corresponding changes in the crystal carrier. The development is 
presented in chronological order from the simple pressure block method to the more elaborate 
forms using automatic truing. The improved equipment is capable of producing 0.001-inch- 
thick quartz crystals and can be used for producing equally thin wafers from a variety of 


other materials. 


I. Introduction 


In order to manufacture erystals that are thin 
enough and have the proper contour to give the 
high fundamental frequencies now being sought, 
it has been necessary to modify the usual lapping 
procedure and to design appropriate apparatus 
Various modifications have been tested; the most 
desirable features have been retained and obvious 
weaknesses corrected (as far as possible) in evolv- 
ing each subsequent model. The apparently de- 
fective units have not been abandoned ; reexamina- 
tions and alterations in the light of added experi- 
ence have been made whenever time and funds 
have permitted. The most successful machines 
develop: din this laboratory produce crystals 0.00] 
in. thiek. Further improvements are expected, 
and the ultimate thickness will depend more on 
problems related to the handling of such fragile 
crystals than on the mechanical limitations of the 
finishing equipment 


II. Conventional Method 


Usually blanks to be lapped are placed in the 
upertures of a metal or plastic nest, which carries 
them through a curved path between the laps 





two heavy flat plates charged with abrasive. A 

modified drill press drives the nest by means of an a 
adjustable eccentric and a bushing lig l Since a = ' “i . 

the nest must be thinner than the crystals to per- 

mit their abrasion, the thickness of the nest is the with such a thickness have maximum fundamental 
limiting factor Beeause of the stretching and frequencies of not more than 14 and 20 me for AT 
buckling that occur, this limit is about 0.005 in and BT cuts, respectively. These cuts are most 
for most available materials, and quartz crystals generally used for high-frequency oscillators 
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III. Cemented Block Method 


The final thickness of the crystal was made inde- 
pendent of the nest thickness by putting the plane 
of pressure application below the upper surface of 
the nest. In the first modification this was ac- 
complished by cementing a *-in.-high by ':-in.- 
square steel block to each ‘:-in.-square crystal 
The cemented units were placed in a nest (0.025 
in. thick) and driven in the usual curved path over 
the lower lapping plate only, the individual blocks 
furnishing pressure previously supplied by the top 
plate Blanks of less than 0.002-in. thickness were 
produced in this Inanner, but they were too wedge- 
shaped to oscillate 

This method of mounting had important dis- 
advantages. Each ervstal had to be mounted 
parallel to the block edges, without strain, precisely 
in the center of the base with the normal to the 
geometric center of the major crystal surface 
passing through the geometric center of the block 


Such conditions were difficult to more than approx 


imate. Furthermore, it was necessary to detach 
the erystal to check thiekness, contour, and fre- 
queney For these reasons other schemes were 
tried 


IV. Free Block Method 


In the next experiment the same apparatus was 
used, but the block was permitted to rest freely 
on the ervstal (fig. 2 Slight movements of the 
block relative to the ery stal constantly redistrib- 
uted the abrasive mixture between the block and 
crystal, thus preventing accumulations that would 
cenuuse unequal distribution of pressure by the block 
However, wedge-shaped crystals again resulted 
ana crystals tended to slip under the nest, so this 
procedure was suspended while other paths were 
pursued 

During the reinvestigation it was observed that 
the movements of both blocks and crystals were 
violent and irregular at the lowest pulley speed, 
so the motor speed Wits reduced by ears of il 
variable autotransformer until, at 40 rpm, the 
Regularly 


at quarter- or half-hour intervals each block was 


motion appeared to be satisfactory 


turned 90° with respect to tts crystal in order to 
compensate for minor variations, such as the over- 
hand of the block. Crystals 0.011 in. in thickness 
were brought down to 0.0025 in. without breakage. 


using the relatively coarse size 400 silicon carbide 
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FIGuri 2 Free blocl appa atlus 


To keep crystals from slipping beneath the thin 
(0.020-in.) zine nest, small weights were placed 
around the loaded openings A heavier nest 
would have made this unnecessary 

Crystals so processed were almost flat. being 
only slightly and equally high at the corners. The 
high corners were removed by edge-grinding to a 
lesser dimension 

Major advantages of the method are: The 
apparatus is simple; the lapping is so gentle that 
breakage is virtually eliminated; individual erys- 
tals may be removed from the nest easily, per- 
mitting more frequent and more satisfactory in- 
spection; crystals of widely varying thicknesses 
may be simultaneously lapped. On the other 
hand, the rate of lapping is slow, and the ery stal 
faces must be parallel at the start, as the block 
is resting on and parallel to the upper surface 
of the ervstal rather than to the lapping plate 


V. Square Block and Cell Method 


As a means of more effectively controlling the 
movement of the crystal and pressure block rela 


tive to each other, both were confined in the -1N) 
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Figure 3 Square block and cel parts 





kieture 4 Square Llack and cell apparat as assembled 
ile the parts shown in the previous figure have been assembled and ar 
t the nest over the i Block n back 


square opening made through the center of a 
l-in. square by ‘4-in.-thick steel plate. This 


unit Was then placed ina nest that had pentagonal 
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openings suitable for lapping inch square blanks 
Counterweights were placed in the other nest 
openings to give more uniform motion 

Because the abrasive that worked into the nar- 
row clearances between the block and plate caused 
binding, the opening in the plate was enlarged 
to about ® in. square. A thin zine sheet with a 

-in. square hole was cemented to its upper 
surface to hold the pressure block in- position 
while permitting the erystal to move relative to 


it (fig. 3 Channels across the block surface 
facilitated this movement by reducing the area 
of adhesion between it and the erystal. This 
modified cell assembly was lapped in the same 
nest and manner used above as shown in figure 4 
Crystals produced by this method were wedge- 
shaped However, the process seemed to hold 
so much promise that a modification that elimi 
nated binding and permitted automatic truimg 
was designed. This unit, which employs a eylin- 
der and plunger, has been put into operation 
recently. Its performance ts described later in 
the article. 


VI. Pentagonal Block and Ring Method 


In a somewhat more elaborate apparatus (fig 
5), pentagonal blocks were rigidly attached to a 
lapped ring. The blocks fitted snugly in’ the 
openings of the nest, which was driven by the 
eccentric as before. Beeause this arrangement 
permitted truing of the pentagonal blocks until 
the lower surfaces were coplanar and parallel to 
the lap, wedge-shaped crystals could be cor- 
reeted to parallelism, an important feature 
Adhesion between blanks and blocks was reduced 
by channeling the under surface of the blocks 
Springs between the ring and nest forced the nest 
into closer contact with the lap, thereby eliminat 
ing the breakage due to erystals slipping beneath 
the nest Weights could be used in place ol 
springs but have not vet been tried 

In this method, as in all the modified procedures, 
the vreatel part of the abrasion took place on the 
under surface of the ery stals, the top surface being 
more polished than abraded. To insure equal 
wear and similar contours on both surfaces, the 
crystals were inverted ut regular intervals Since 
the contours of the cry stals were dependent on the 
trueness of the pentagonal blocks, the blocks (still 
attached to the ri w) were lapped before and at 


suitable intervals during each run. 
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Fieure 5 Pentaqonal block apparatus 


Res ilts achie ved by the use of the pentagonal 


block svstem were encouraging Deviations from 


flatness were radial rather than wedge-like Th 
variation did not exceed 0.00004 In. and was. for 


the most part concentrated ina rim-like burden 


the enelosed area being very nearly flat and 


parallel BT ervstals with a thickness of about 


O.0025 and oa fundamental frequeney of 45° me 


Waele produced Sueh rvstals could bye etched 


toeven th ohet req rencies 


Experience with the pentagonal block svstem 


suggested two modified procedures In one the 
nest was eliminated mn the other uv thick nest 
vith « tlar hol Wis used 


VII. Driven Ring Method 
Phe re 


ipon which were mounted 


st-less upparatues consisted of a flat ring 
blo ks 


each with an enetrelit ring or collar to retain 


eviindrical 
the ervstal (fig. 6 The collars and the under 
channeled to improve 


surfaces of the blocks wer 


ireulation of the abrasive mixture An eccen- 
‘ netting through foul rigid spokes, drove this 


assembly over the lower lap Because the ring 
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ng apparatus 
,;api 


and spokes were of thick steel, the load on the 
crystals was quite heavy. As a result, crystals 
began to break when approximately 0.004 in 
thick. The investigation has been discontinued 


until a ring of lighter metal can be made 


VIII. Ring and Thick Nest Method 


This modification made use of the same type of 
ring as above but without the spokes, since thy 
ring was carried by the nest (fig. 7.) Cylindrical 
plugs were used without collars. Instead, holes 
were carefully drilled and reamed in a thick and 
rigid nest. Each hole was a sliding fit on any 
plug of the ring Indeed the main advantag 


of using the circular rather than the pentagonal 
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openings lay in the ease with which apertures 
and plugs could be made to close tolerances 
As in the pentagonal block method, the nest 
served both to retain the erystals and to adjust 
the pressure on them. This latter effect was 
accomplished by the use of small conical springs, 
the apex ends of which fitted tightly over studs 
located midway between the lapping blocks and 
the base ends of which pressed the nest into 
closer contact with the lapping plate. Changing 
the number of springs adjusted the pressure on the 
erystals, permitting heavier pressure and faster 
lapping on the thicker crystals and similarly, 
lower pressure and slower lapping on the thinner 
ones. Cross-channeling the under side of the 
nest with grooves ‘45. In. wide by 6 In. deep at 
in. intervals allowed freer circulation of the 
abrasive and reduced adhesion between the lap 
and the nest. 

Although this method is still under investiga- 
tion, some conclusions have been drawn. Since 
the nests were thicker than those available com- 
mercially, the apertures maintained their size and 
shape much better than did the openings in the 
thin nests. Consequently, the clearance between 
each plug and its respective hole remained small 
and the resulting rim was much less pronounced 
In the case of square crystals, the rimming was 
confined to the corners where its effect on per- 
formance appeared to be negligible. So far, 
crystals have been ground to a thickness of 0.002 
in. by this method. 

Although no definite limit was imposed on the 
final thickness of erystals produced with the 
equipment last described, the hazards connected 
with removing the ring and nest in order to inspect 
crystals less than 0.002 in. thick made another 
method seem more desirable for the final stage of 


lapping 


IX. ‘Inkwell’’ Method 


As mentioned earlier, a cylinder and plunger 
modification has been designed on the basis of 
experience with the square block and cell unit. In 
the improved form (fig. 8) the retaining cell con- 
sists of a cireular aperture in a heavy-walled 
brass block, approximating an inkwell in size and 
shape. The conical shape is a compromise be- 
tween height and stability, and the base is broad 
to maintain its true surface. The bottom part of 
the block is a detachable annular disk, the central 
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FIGURE 8. Inkwell apparatus parts, 


he keyed plunger, conical body, and annular disk are shown assembled 
t the left. In the foreground may be seen the upper surface, chamfer, and 
bore of the disk part Phe lower surface and its relation to the bore ‘ lent 
n the unit slight to the right, while the plunger may be seen in the right 
backg ind 


hole of which contains the erystal during lapping. 
The diameter of this hole is just enough larger 
than the plunger to prevent binding by the abra- 
sive. For the same reason the cell is quite shallow. 
A recess on the upper block together with the 
chamfered upper surface of the disk prevents 
abrasive, from reaching the more closely machined 
upper bore. Since the disk is detachable, it can 
be replaced easily when wear renders it unusable 

The evlinder bore is a smooth fit for *-in. drill 
rod. For lighter pressure, hollow) aluminum 
plungers are being used. Close fit) (especially in 
the upper bore) is essential to keep the plunger 
from wobbling. The plunger and block base are 
trued with respect to each other against the lap 
A pin in the block wall engages a keyway in the 
upper part of the plunger thereby preventing tts 
rotation while permitting vertical motion. In this 
manner the trued surfaces of the plunger and base 
remain parallel except for the factor of prolonged 
wear). [It is customary to true these surfaces at 
about S-hour intervals. As in previous lapping 
units, the working surfaces are channeled to 
improve circulation of the abrasive and to reduce 
adhesion. 

To operate, the loaded unit) (cell block with 
plunger and erystal) is placed in one of the circular 
apertures of a heavy nest and carried over the 
lower lap in the usual curved path (fig. 9). Since 
the aperture is 0.015 in. larger than the cell block 
and the cell aperture is considerably larger than 
the enclosed ery stal, all three cry stal, block, and 
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hest execute different eccentric motions The 


present lap and nest can accommodate eight units 
The erystal in any one of these units can be easily 
examined without disturbing any other crystal, a 
verv useful feature 

[ sing this apparatus « rvstals have been lapped 
to about 0.001 in. with practically no breakage 
The edges were smooth and free from = radial 
eracks, and the faces were, in’ general, satis- 
fae torily flat Naturally 
fragile and handling them in the usual routines of 
measurement, plating, and 


cleaning mounting 


These 


problems, rather than the method of lapping 


gives rise to important new problems 


may set the ultimate limit for ervstal thickness 


X. Tall Plunger Method 


With the intention of regulating the piston 
wobble, a modification of the above unit has just 
been completed. In the new unit the plunger is 
taller and bears on rounded points rather than on 
a wall. Twelve screws with rounded ends serve 


as the bearimg points, six being equally spaced 


nround the piston hear its top anda anothe SIN 


erystals so thin are 


FIGURE 


Afte I 


these screws have been adjusted to permit sliding 


similarly spaced near the bottom (fig. 10 


without lateral play, the transverse serews ar‘ 
tightened to maintain the fit. This modification 


has vet to be tested. 


XI. Conclusion 


Of the ten methods deseribed the last) thre 
Although they have 


been developed primarily for the production of 


were the most satisfactory 


thin quartz oscillator plates they can be used for 
lapping a variety of mineral and ceramic materials 
For example, plates of a ceramic dielectric have 
been lapped by the inkwell apparatus to 0.0015-in 
thickness. In such experimental work the indi 
vidual plunger and block units are especially 
useful, since only a minimum of the test material 
is needed. As previously mentioned, a numbet 
of such units may be used with one driving nest 
various thicknesses 


to lap test specimens of 


simultaneously Hlowever, since the driven nest 


with ring-mounted plugs accommodates more 


ervstals on the same size lap it is more economical! 
When many plates of the same thickness are to b 
made. Both methods have their special features 
and it is likely that development will continu 
along two paths to provide machines that will 
better serve, on the one hand, for production of 
experimental or special orders, and on the othe 


hand, for commercial quantities 


W \sHINGTON, June 14, 1949 
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Instability in Shear of Simply Supported Square Plates 
With Reinforced Hole 


By Wilhelmina D. Kroll 


The critical shear stresses were determined for five elastic, simply supported square 


plates with central circular holes of diameter one-eight and one-fourth the length of a side. 


Comparison of the numerical results obtained with those computed for plates without holes 


shows that, although an unreinforced hole may 


cause a large reduction in the critical shear 


=tress of the plate, reinforcement of the hole by a circular doubler plate causes a substantial 


nerease in the shear buckling load 
I. Introduction 


Holes must frequently be cut in the stressed 
skin surfaces of airplane wings or fuselages to 
provide access to the interior of the wing or fuse- 
lage. To prevent weakening of the entire strue- 
ture thereby, the hole is usually reinforced by a 
circular doubler plate riveted to the inner surface 
of the sheet 

The present report is one of a series of studies 
Navy 


Department, in order to understand more fully 


undertaken for the Bureau of Aeronautics, 


the effeet of reinforcements in strengthening and 


stabilizing the structure around a hole. In 


reference 1} a theoretical Investigation Was 
presented of the effect of various reinforcements 
on the stresses and displacements near a small 
circular hole in a plane sheet under uniform ten- 
sion in all directions. Numerical values of stress 
concentration and of radial displacement at the 
edge of the hole for various doubler plates indi- 
cated that the greatest reduction in stress con- 
centration for ai given volume of reinforcing 
material occurs when the outer radius of the 
reinforcement is smallest and that the reduction 
in distortion of the hole depends principally upon 
the volume of the reinforcement 

made on plates with reinforced 


Tests were 


circular holes, under uniform tension ino one 


direction only, reference [1], to check the theoreti- 


cal analysis, and it was found that, if the reinforce- 
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ments were attached to the sheet with two con- 
centric rows of rivets or bonded to the sheet, the 
theoretical analysis gave a sufficiently accurate 
description of the stresses and displacements in 
the neighborhood of the hole. 

In reference |2], results of compressive tests of 
curved panels with circular holes unreinforeed or 
reinforced with doubler plates are given, which 
show that curvature has no appreciable effect on 
the median fiber strain distribution around a hole 
and that the plane stress theory can be used to 
reliably predict the stress distribution in the un- 
reinforced regions of curved panels with holes. 

A theoretical analysis for the stress distribution 
im a flat 
loaded by a pin was presented in reference [5]. 


plate near a reinforced circular hole 


Results obtained in the theoretical analysis and 
those from a test of a plate of sandwich con- 
struction showed good agreement. 

A method is presented in reference [4] for deter- 
mining the instability under edge compression of 
a rectangular plate with a reinforced circular hole 
The values of critical compressive stresses of 
square plates obtained by this method indicated 
that, although the buckling stress of a plate is re- 
duced only a small amount by the presence of an 
unreinforeed hole, reinforcement of the hole causes 
a substantial increase in’ the buckling. stress 
Reinforcement of a hole, therefore, can be expected 
to stabilize as well as strengthen the structure in 
the neighborhood of a hole when the structure ts 


under compressive load 














Using the method of reference |4], a study is 
made in this report of the effect of the reinforce- 
ment of a circular hole in a square plate on the 
stability of the plate under shearing loads uni- 


formly distributed along the edges of the plate. 


II. Method of Analysis 


In reference [4], an energy method was presented 
for computing the compressive buckling load of a 
simply supported elastic rectangular plate having 
a central circular hole reinforced by a circular 
doubler plate. This method can also be used to 
determine the critical shearing stress for simply 
supported rectangular plates having a central 
circular hole 

The integrals 7; and /, ,;which must be evaluated 


for the plate are 


b= | [Doe ty) 
surface 


ow ow o-w \~ 1) 
2 l Md —= d. dr 
S. OY (ss) | es 


=f Jo[% (2) 4 


surface 


0 ow 
(ee 


where hh 





T ow ow 
) TZ da dy 
T. OF OY . J 


plate thickness, function of « and y; 


rectangular coordinates with origin at 


r. 

center of plate; 

w—lateral deflection of plate; 

D= EN 12(1-n*), flexural rigidity of plate, 
function of « and y: 

u- Poisson’s ratio, 0.3; 

o,— tensile stress in « direction; 

a, — tensile stress in y direction; 


stress ratios 


shear stress in «and y directions: 


shear stress in s and y directions far 


from the hole 


G./\Te, GuitTe and Txy Ta 


just 


prio! to buckling may be obtained from reference 
4] in the following w TA If the origin of coordi- 
nates is taken as shown in figure 1, shear stresses 
along the edges of the plate will be equivalent to 
tensile stresses acting at the edge of the plate ut 
#—45° and equal compressive stresses acting at 
the edge of the plate at @ 15 
F and A of reference 1] are 


The constants 


zero in the ease of 
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shear loading. stresses ¢,/Te, O@/Ta, and 


The 
rT referred to polar coordinates r, @ are then 


riveh as 


o,/T A+3CR'/r'+-2DR?/r 7 
cos 2 (64 1) 
O9/T« A+6Br P+ 3CR 
> 2) 


cos 2 (64 4 
T109/T A+3Br/R 


3CR'/r' 
DR? /r*) sin 2( 64 1) 





where 7, @=polar coordinates with origin at 
center of hole (see fig. 1 
#—angle between radius r to point and 
positive horizontal axis. Positive 
#@ measured in counter clockwise 
direction; 
R—radius of hole; 
A, B,C, D coefficients with different values in 


sheet and reinforced region 














\ # 
N A 


a 











Fiagure 1. Subdi nto three portior 1, B 


and C, in evaluating 1, 


ision of plate 


The values of A, B, C, 
as outlined in reference |4] The stresses in rec- 


and i) nia be det rmined 


tangular coordinates may be computed from the 
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Stresses in polar coordinates, eq 2, by use of the 


conversion formulas: 
a a, Cos*6+- a9 sin® 6—r~ sin 26 


a, sin"6+ a, cos* 6+ 7,9 sin 20 
Oe 
r, sin 26+rw cos 26 


The lateral deflection, w, of a simply supported 
rectangular plate of length @ and width 6 can be 
approximated by the first terms in the trigono- 
metric series 


») » 


Td ry - ona . 2a 
w=), COS COS 7 doo sin sin --+t 

a 4 r a b 
wr Ty Omer Ty 

COs cos “ = fle, COS cos ; + 
a b : a b 

. 2ar . 4ay 

(og SIN sin ere 4) 
bh 


After substituting eq 4 into eq 1 and performing 
the indicated integration by a numerical proce- 
dure, the integrals /; and /, become quadratic 
expressions in terms of the coefficients a,,, dos, 
Qj, . . . . The eritieal value of r., the 
shear stress far from the hole, at which buckling 
of the plate occurs is that value of +. which re- 
duces to zero the determinant of the coefficients 
Of @), dy, Ay, ete. in the set of simultaneous 


equations 





ol ol; ) 
-T. 0 
On; Od); 
° a) 
Od, o Ol os 5 
Oa Od, 
ee ae a ee J 


Ill. Numerical Integration 


The evaluation of the integrals in eq 1 over the 
surface of the plate could not be done directly, 
since the plates have a rectangular outer boundary 
and a circular inner boundary and involve a stress 
that is a complicated function when expressed in 
rectangular coordinates. 

The integral J,, eq 1, was taken as 


Nh=Tattho—he (6) 
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where /,, =the integral /, for the surface enclosed 
by the outer rectangular boundary 
of the plate, taking the flexural 
rigidity D=D,, the value of Din 
the unreinforced portion of the 
plate; 

/,,—the integral /, for the surface enclosed 
by the outer circular boundary 
of the reinforcement — taking 
D—D,—D,, where D, is the value 
of )) in the reinforcement; 

/, the integral /, for the surface enclosed 
by the circular boundary of the 
hole taking dD D),. 


The double integration necessary to evaluate 
I, could always be done directly, but this was 
not possible in evaluating J,, and J,,. The inte- 
grals J,, and J,;, were obtained by integrating 
directly with respeet to and using Gauss’ method 
of numerical integration, reference [5], for inte- 
grating with respect to y. In all cases except 
case 2, five Gauss points were used in this numer- 
ical integration as computations showed that, 
in inereasing the number of points used from 
three to five, the value of the critical stress was 
changed only 2 percent. If more than five points 
had been used, the change in the value of the 
critical stress would have been even less. For 
case 2, three Gauss points were used in evaluating 
T,. beeause the circular area over which /,. was 
evaluated had a much smaller radius for case 2 
than for the other cases 

The integral J, was evaluated as the sum of 


three integrals 


To= Toqt- Lon t Lec, 7 
where J, 


the integral /, for the cireular disk 
of the reinforced area (A, fig. 1 

/,,=the integral 7, for the circular disk 
between the outer boundary of 
the reinforcement and the largest 
inscribed cirele in the plate (B, 
fig. 1; 

I, the integral I, for the remainder of 


the plate (C, fig. 1 


The integral for each circular portion was deter- 


mined, using Gauss’ method of numerical inte- 
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gration, by first integrating numerically in a 
circumferential direction and then in a= radial 


direction 


plate Was obtained also by using CGiauss’ method, 


first integrating numerically in the z-direction and 


then in the y-direction 


| kdA e [(0.0481444F;, + 0.07703 LOF, + 0.0481444 F, + 0.0380048 F + 0.057 LS68F;, 4 


0.0914988F), + 0.057 L868 F; 4 


0.0337 158F), + 0.0539452F 
0.0142967 F O.OLSIS805F)-4 


0.0134864F,, + 0.00842899 F,,] 





The integral for the remainder of the 


0.0727220F, 4 
0.0337 158F; 


0.0290888 F). 4 
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A typical distribution of Gauss points for com 
puting /, is shown in figure 2 for one-eighth of a 
square plate with a central circular hole. For 
this distribution, the integral of a function F over 
the portion of the plate surface enclosed by heavy 


lines in figure 2 is 


0.116355F,+-0.0727220F,, 


+ 0.0142967 F\, + 0.02287 47 F), 


O0.O1ISTS05F), + 0.00842899 F,, 





where F), fo, Fy... . Fao are values of the function 
Fat the positions shown in Figure 2 

In each case a sufficient number of Gauss points 
was used to reduce the estimated error to less than 
5 percent. Twenty-two points for each octant of 
the plate were used for cases 1, 3, and 4; twenty- 
eight points for case 2; and twenty-five points for 
cuse 5 

As an indication of the adequacy of the nume: 
ical integration methods used, the critical shea: 
stress of a square plate without a hole was deter- 


mined by approximating the deflection by 


Ta Ti 2rs Jry 
uw a cos oon -" +6 sill sin ’ 
‘ a i a a 


integrating for /,;, eq 1, exactly and integrating 
for J,, eq 1, both exactly and by numerical 
integration with 22 points in one-eighth of the 
plate as shown in figure 3. The resulting critical 
stresses differed by O.S percent In the case of a 
plate with a reinforced hole, it is probable that 
the more complicated stress distribution and th 
use and prominence of higher order terms in 
the series used for the deflection will cause thi 


error to be somewhat higher. 


IV. Convergence of Deflection Function 


The correct value of +. for buekline of th 


plate is that value of 7. which reduces thy 
determinant of the coefficients of a),, do, ay3, 

in the infinite set of eq 5 to zero. In order t 
limit the work of computing 7. to a finite amoun 


preliminary computations were made to dete! 
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were most important in evaluating r+. 


3. table 1. 
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/ ‘ case ol plate vithout a hole 


Using the distribution of points shown in figure 


1 to carry out the numerical integration for J,, 
was IS.7S teh a 


only the first two terms of eq 4, 


the eritieal shear stress + using 


Td ry 2 ra ry 
COs cos ra cos COS . 
a a s a a 
as the deflection function 
Td TY Jr Jr 
dy, COS cos ~T da. sin sinh T 
( a i 
Td sry sr TH 
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mine which terms in the deflection function, eq 4, 
The 


computations were made for the plate of case 
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Figure 4. Dis 


bution of Gauss points for 


dy of convergence of de flection f 
T was 18.10 /h? a*, a deerease of 
Using 
- ) 
Td wy = Wd 
we a Cos COS Ta siti 
( a ad 
Ord ry 
ad COS COs » 
a a 
T was 17.46 Teh ay, a decrease of 
Using 
Td a] 2rd 
u a COs ‘oa — a Sih 
a a 
Jr Ty 
(log Sh sin “ tTdg Sil 
ad * 
r. Was 17.37 Fh? a?, a decrease of 
Using 
Tl Ti Jari 
u a COs cos a sill 
a a 
Td ries aye 
a cos cos ~ ———_ COS 
‘ a a 
was 18.50 Kh? a a decrease of 
Using 
Td wij Dr. 
u ad COs COs . a sil 
‘ ad ad 
Td 47r1/ 
(4, Sin Sih ’ 
a a 


was IS.5S Ih? a 


_ = deerease ol 


| 


‘ valuat 


snection, 


5a 


ng Isin 


3.0 pereent 


. 2ry 
sill 
ad 
7.0 per 
2ry 
sin 
ri/ 
Sin 
ad 
4.0 pere 
O. 
sin 
ad 
rh 
COS . 
] 5 per 
— 
wh 
Sin 


lL percent 


+ 


‘ent 


ent 


469 








The use of additional terms of the trigono- 


metric series with the first two terms probably 


would have shown a proportionately smaller 


decrease in the buekling load. Therefore, it is 


believed that the function 


Td ri Jr 2ry 
u a COS cos * +a sinh sin T 

a a a 
Td Ti J ta 

(lag COS COs “ = dog SIN Sl T 
a a a 
tors 2ry 

dg sin ’ 
( a 


approximated the lateral deflection of a square 
plate with simply supported edges under shearing 
loads with sufficient accuracy. This deflection 
function was used in the computations of the 
critical shearing stresses 

After the eritieal shearing stresses had been 
computed for the five square plates considered in 
this report, the relative importance of the various 
terms in the deflection function was determined. 
Values of a. coetfi- 


CIENUS Cga, Coy, and dy 


were 25 to 42 percent of ay: 
were 0.4 to 5.0 percent of 
a It seems reasonable to assume that higher 
order terms would have been even smaller as 
compared to @, and that the most important 
terms were used for the deflection in all cases. 

The number of points used in the numerical 
integration for 7, has an effect on the accuracy of 
the determination of the critical stress. To 
earry outa complete Investigation of convergence, 
as more points are used, was not possible because 
of the amount of work involved. To obtain a 
partial indication of the accuracy to be expected, 
however, the critical stress for case 3, table 1, was 
computed using various deflection functions and 
using both the distribution of points in figure 2 
and that in figure 4 for computing /,. The results 


are given in table 2. 


Tapie 2. Values of critical stress obtained from the spacings 
of Grau points cnvestiqated 
( oe ( r 
leet 1 I Ditter 
i4 figure 4 M rl I 
Pe 
STS ED a] hea 
is BRD a 4.56 EA2/a ‘ 
EA S21 FA 
S50 Eht/a 58 Eh 6 
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It was decided on the basis of these computa- 
tions that the use of the point distribution in 
figure 2 would be acceptable. 

The over-all accuracy of the computations of 
critical stress are estimated to be such that the 
critical stresses are within 7 percent of the exact 


values for the cases investigated. 
V. Results and Discussion 


The critical shear stress was determined for five 
square plates with reinforced and unreinforced 
holes, figure 5. The dimensions for the plates are 
given in table | 

The analysis gave the critical shear stress fat 
from the hole r.., corresponding to the stress dis- 
tribution for a plate under uniform shear at an 
infinite distance from the hole. The shear stress 
tr» IS larger than the average shear stress on the 
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FIGURE 5 Loading of square plates with s mply support 


edges (Dimensions in table 1 
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center cross section obtained by dividing the load 
on that section by ah. A better value of the 
average buckling stress 7,, of the plate was ob- 
tained from 

1 P+ P, Q) 
ad 2 ah ; 
where P, shearing force acting on edge of plate; 
P,—shearing force transmitted across center 

section, y= 0. 


The value of P, was obtained by numerical integra- 
tion and the value of P, by direct integration of 
the stresses r,, given by eq 3. Values of r., 


T P, ah, T2 


P,/ah, and +r,, are given in table 3. 


TasLe 3. Critical shear stress for plates investigated 


C use T T Tt T 
rh 
5.46 — 
Eh Eh Eh _ Bh 
6.0 out 6.11 O58 
: a a 1 
Eh Eh Eh Eh 
y 4 s s es] - sw 
1 a? 
Eh _.. Eh _ Eh _ Bh 
4 | re 17.01 2% 
t a ! 
Eh Eh ., Eh Beh 
i Wd 10.46 rh.) W12 
! a t 
Eh Eh . Fh Eh 
A Jas 27 wD a 
t a ! 


It was assumed that the boundaries of the 
square plates were sufficiently far from the hole 
so that only a uniform shear was acting on the 
boundaries. Numerical integration of the stresses 
o, as given by eq 3 on the boundary of the plate 
showed that there was an average compressive 
stress on half the side of the plate and a tensile 
stress on the other half equal to 3 to 12 percent 
of the shear stress. This accounts for the small 
differences in the values of 7; and 7. 

The shear buckling stress of a simply supported 
square plate of constant thickness is given on 
page 360 of reference [6] (taking Poisson's ratio 


us 0.3) as 


Ta=EA6 kh? a’. (10 


Comparison of this value, case 0, table 3, with the 


values of the critical stresses for cases 1 and 2 


Instability of Plate with Circular Hole 


shows that the unreinforced holes of diameter one- 
fourth and one-eighth the length of a side reduced 
the buckling stresses by 22.6 and 0.2 percent, re- 
spectively. Comparison of the value for case 0 
with the values for cases 3 to 5 of table 3 shows 
that reinforcement of the hole increased the shear 
buckling stress of the plate over that for a plate 
without a hole from 20 to 236 percent. These 
percentages are to be taken only as indications of 
the order of magnitude since the accuracy of the 
present determination of critical stresses is esti- 
mated to be 7 percent. 

The effect of thickness of the reinforcement on 
the buckling load is seen by comparing cases 3 
and 4, table 3. In these cases, the hole size and 
area of the reinforcing material is the same, but 
the thickness of the reinforcement of case 4 is only 
one-half the thickness of the reinforcement of 
case 3. The buckling load for case 8 is 71 percent 
higher than that for case 4. 

The effect of shape of reinforcement on the 
critical shear stress is indicated by comparing 
The volume of material 


cases 3 and 5, table 3. 
reinforcing the holes in these two cases is the 
same, but the reinforcing material is concentrated 
nearer the edge of the hole in case 5 than it is in 
case 3. The buckling stress is 64 percent greater 
for case 5 than for case 3. 

Comparison of the values for case 1 and cases 
3 to 5 shows that reinforcement of the hole raises 
the critical shear stress from 55 to 334 percent, 
depending on the thickness and shape of the 
reinforcement. 


VI. Conclusions 


Using a numerical procedure for evaluating the 
integrals for the energy stored in a plate, the 
critical shear stress can be estimated for plates 
with circular holes and doubler plate reinforcement. 

The critical shearing stress of plates may be 
reduced considerably by the presence of unrein- 
forced holes. A square plate with a hole diameter 
one-eighth the length of a side showed practically 
no reduction in critical shearing stress, but a square 
plate with a hole diameter twice as long showed a 
reduction of 23 percent. 
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Reinforcement of the doubler plate type causes 
marked increases in the critical shear stress. The 
critical shearing stresses were 55 to 334 percent 
higher than that for a plate with the same size 
unreinforeed hole 

The volume and shape of the reinforcement has 
a marked effect on the critical shearing stress 
The critical stresses for two plates with the same 
size hole and same area of reinforcement, but 
different thickness, were compared. The critical 
stress increased from 10.12 /h?/a? to 17.27 Eh?/a?, 
an increase of 71 percent, when the thickness of 
the reinforced region was increased from 1.5 h to 
20h 
the same size hole and same volume of reinforcing 
The plate with the 


The critical stresses for two plates with 
material were compared 


material concentrated closer to the edge of the 


hole had a 64 percent larger critical shearing stress 
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Physical Properties of Sixty API-NBS Hydrocarbons 


By Alphonse F. Forziati * and Frederick D. Rossini 


The results of experimental determinations of the density at 20°, 25°, and 30° C, the 


refractive index, np, at 20°, 25°, and 30° C, 


and the boiling point and pressure coefficient of 


the boiling point, at 1 atmosphere, are reported for 60 hydrocarbons of the API-NBS series, 


including 17 paraffins, 14 alkvlevclopentanes, 8 alkvleyclohexanes, and 21 alkylbenzenes. 


I. Introduction 


As part of the work of the American Petroleum 
Institute Research Project 6 on the fractionation, 
analysis, purification, and properties of hydro- 
carbons, measurements are being made of the 
physical properties of the API-NBS_ series of 
highly purified hydrocarbons. This paper gives 
the results of measurements of refractive index, 
density, and boiling point of 60 API-NBS hydro- 
carbons, including 17 paraffins, 14 alkylevelopen- 
tanes, S alkvlevclohexanes, and 21 alkvlbenzenes 


II. Compounds Investigated 


The compounds whose properties were measured 
in the present investigation were samples from 
the API-NBS series of highly purified hvdro- 
carbons, which are being prepared through a 
cooperative undertaking of the American Petro- 
leum == Institute and the National Bureau of 
Standards 

These samples of API-NBS hydrocarbons have 
been made available by the American Petroleum 
Institute and the National Bureau of Standards 
through the APL Research Project: 44 on the 
“Collection, Analysis, Calculation, and Compila- 
tion of Data on the Properties of Hydrocarbons” 
The samples were purified at the National Bureau 
of Standards by the APL Research Project 6 on 
the “Analysis, Purification, and Properties of 
Hvdrocarbons”, from material supplied by the 


following laboratories: 


1} t tion w performed at the National Bureau of Standard 

rt of tl k of the American Petroleum Institute Research Project ¢ 
the Arma . Purification, and Properties of Hydrocarbons 

Res hA ciate on the American Petroleum Institute Research Pro 


| Bureau of Standard 


Physical Properties 


n-Heptane, 2,4-dimethylpentane, —-nonane, 
2,2,5-trimethylhexane, benzene, 1,2,3-trimethyl- 
benzene, and 1,2,4-trimethylbenzene, by the 
American Petroleum Institute Research Project 
6 at the National Bureau of Standards. 

2-Methylhexane and 3-methylhexane, by the 
Ethyl Corporation, Detroit, Mich. 

3-Ethylpentane, 3,3-dimethylpentane, 3,3-ci- 
ethylpentane, 2,2,3,3-tetramethylpentane, ¢/s-1,3- 
dimethylevelopentane, — trans-1,5-dimethyleyelo- 
pentane, n-propyleyclopentane, —isopropyleyclo- 
pentane, 1,1-dimethyleyclohexane, -propylevelo- 
hexane, isopropyleyclohexane, 1,1,5-trimethyley- 
clohexane, — v-butyleyelohexane, — isobutyleyelo- 
hexane, sec-butyleyclohexane, tert-butvleyvclohex- 
ane, 1,3-dimethvlbenzene,  1,4-dimethylbenzene, 
l-methyl-2-ethvlbenzene, l-methyl-3-ethylben- 
zene, 1,3,5-trimethylbenzene, v-butylbenzene, iso- 
butvibenzene, sec-butylbenzene, and tert-butvl- 
benzene, by the American Petroleum Institute 
Research Project 45, at the Ohio State University, 
Columbus, Ohio, under the supervision of C. E. 
Boord. 

2,2-Dimethylpentane, by the Socony-Vacuum 
Laboratories, Paulsboro, N. J. 

2.3-Dimethyvlpentane and 2.2.5-trimethyl- 
butane, by General Motors Corporation, Detroit, 
Mich 

2.4.4-Trimethylhexane, — 1,1-dimethyvleyelopen- 
tane,  ¢/s-1,.2-dimethvleveclopentane, — frans-1,2- 
dimethyvleyelopentane, —1,1,2-trimethyleyclopen- 
tane, cis,cis trans-1.2,4-trimethyvlevclopentane, 
cis transcis-1,2.4-trimethvleyvclopentane, —l-meth- 
vi-l-ethvleyvclopentane, and l-methyl-e/s-2-ethyl- 
evelopentane, by the Hydrocarbon Laboratory, 
Pennsvivania State College, State College, Pa 
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2,2,3,4-Tetramethylpentane, 2,2,4,4-tetrameth- 
vipentane, 2,3,3,4-tetramethylpentane, 1-methyl- 
t-ethylbenzene, —1,2-iethylbenzene, and — 1,4- 
diethylbenzene, by the National Advisory Com- 
Flight Propulsion Re- 


search Laboratory, Cleveland, Ohio. 


mittee for Aeronautics, 
1.1,3-Trimethyleyclopentane, by the American 
Petroleum Institute Research Project 45, at the 
Ohio State University, Columbus, Ohio, and the 
Gulf Oil Company Fellowship, at the Mellon In- 
stitute of Industrial Research, Pittsburgh, Pa. 
Ethylevclopentane, by the American Petroleum 
Institute Research Project 45, Ohio State Univer- 
Ohio, and the Hydrocarbon 
Pennsylvania State College, State 


sity, Columbus, 

Laboratory, 

College, Pa. 
Methylbenzene, by the Humble Oil and Refin- 

ing Company, Houston, Tex. 

by the 

Monsanto Chemical Company, Dayton, Ohio 

the Standard Oil 

Development Company, Elizabeth, N. J. 


Ethylbenzene and isopropylbenzene, 


1.2-Dimethylbenzene, by 


n-Propylbenzene, by the Dow Chemical Com- 
pany, Midland, Mich 

1 .3-Diethylbenzene, by the National Advisory 
Flight 


Research Laboratory, Cleveland, Ohio, and the 


Committee for Aeronautics, Propulsion 
American Petroleum Institute Research Project 
$5 at the Ohio State University, Columbus, Ohio. 


The purification and determination of purity 


and freezing points of these compounds are 
described in references |1, 2, 3, 4, 5] a 

Figures in bracket ndicate the terature references at the end of th 
paper 
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It is believed that in each case the impurity was 
of such nature and present in such small amount 
that the properties measured were not affected 
beyond the indicated limits of uncertainty 


III. Refractive Indices 


The measurements of refractive index were 
made with a precision refractometer,’ by reference 
to the three hydrocarbons, 2,2,4-trimethylpen- 
methyleyelohexane, and 


(toluene), for which the values of refractive index 


tane, methyl benzene 
were determined on a spectrometer with hollow 
prism by L. W. Tilton, of the Optical Instruments 
Section of the National Bureau of Standards 
These three hydrocarbons are NBS Standard 
Samples of hydrocarbons certified with respect to 
20°, 25°, and 30° CC. The 
values of refractive index for the three reference 


refractive index at 


hydrocarbons were reported with an accuracy of 
+ 0.00001, and 
methylpentane for the 


were used as follows: 2 2 4-tri- 


paraffin hydrocarbons, 
methylevclohexane for the cvcloparaffin hvdro- 
carbons, and methylbenzene for the alkylbenzene 
hydrocarbons. By this procedure, the refractom- 
eter was actually used only to determine the differ- 
ence in refractive index, np, between the refer- 
compound and the undet 


ence hydrocarbon 


investigation. The experimental results on refrac- 
tive index are given in table 1. Individual measure- 
t 0.00002 to 


The accuracy of the tabulated values, 


ments were reproducible within 
+ 0.00003. 
including the effect of impurities, is estimated to 
be +0.00005 to +0.00008. 


* Modified Abbe-type, Bausch and Lomb 
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TABLE I. Values of refractive index, de nsity, and howling point for sitty API-NBS hydrocarbons 
Pressure Refractive index * up rempera Density rempera 
Boiling coefficient ture coeffi 
. point at of boiling cient of re ture coeffi 
Compound Formula 49 mm__ point at 760 fractive in oie gs 
e mm He 20° ¢ 26° C wc dex (up 2° C 25° C wec ae x “ 
dtdP ut 36° C . 
PARAFFINS 
( Cimmiig Cc g/ml g/ml g/ml g/ml °¢ 
Heptane C-H WS. 427 0. 04481 1. 387M 1. 38511 1. 3825s 0 000506 0.68376 0. 87951 0 0 000850 
2-Methylhexane C-Hy oO. O52 o44sl 1. SS485 1. 38227 1. 37969 OOOSL6 H7S50 67439 H7001 OOO858 
3- Methylhexane CH 1. S50 04459 «L388 OL S809) «1 SN OOO510 OSTIS HN205 67852 OO0N6 I 
}-Ethylpentane CHy 92. 475 04482 «1.39830 1.39084) 1. 3RRDO OOOS10 HONT6 HYSYS CSU4S OO0868 
2,2- Dimethylpentane C Hy 79. 197 (4394 «1.38215 LL. S7955 1. 37605 QuO520 O73S85 HiGAS H0.508 OOONTT 
2,3-Dimethylpentane C-Hy SU. 7M4 (4482 «1.39106 «1.38045 1. S884 OuOs02 64508 H9091 H8050 OOOSSS 
2,4-Dimethyl pentane CrHy SO. 500 04376 «1.38145 1.37882) 1.37619 OOO5 26 67270 HORS2 HHSS3 OOOSST 
},3- Dimethylpentane CH SH. O4 04509 «1.39092 «1.38842. 1. SR5O2 000500 69327 HSGOS OSs OOOH 
2,2,3-Trimethyl butane C-Hy SO) SAD 04484 (1380441 SHO? 38440 000504 6voll HSS8S isi OOOSST 
Nonane CyHy 15). 70S 4967 1.40542 1.40311 1. 40080 ooode2 71708 T1381 70004 QO 76S 
2,2,5-Trimet hylhexane CyHy, 124. 084 OASIS 1.39972 «1. B9T2ZN 1. BO484 OOOASS 70721 TOS22 69905 OOOST6 
2.4,.4-Trimethylhexane CsH, 130. 648 04960 «1.40745 1.40515 1. 40285 000460 72381 T2007 71612 oooTAa 
},3- Dieth ylpentane CyH: 146, 16s O5109 «1.42051 1.41837 1.41623 0042s 5 744 wT 2 
2,2,3,3-Tetramethylpentane Cyl, 140. 274 O5124 1.42360 1.42140 1.41920 Ooo440 TH 74925 OOOT4 
2,2,3,4-Tetramethylpentane (gH, 133. 016 O5030 1.41472 1.41246 1.41020 O00452 TSSOS 73144 OOOTSI 
2,2,4,4-Tetramet hylpentane CyH, 122. 284 (4937 «1.40604 «1.404591. 400224 OOO TO 71947 71S 71167 OOOTSO 
2,3,.3,4-Tetramethylpentane CyH 141. 551 4225 1. 42222 1. 420038 1. 41784 000458 TATS TALS 74740 (OT AS 
ALKYLCYCLOPENTANES 
Ethyleylopentane C-Hy 103. 466 0.04623 1.41981 1.41730 1.41479 0.000502 O. 7047 0 0. 75780 0. OOOS6T 
1,1-Dimethyeyclopentane CoHy 87. M46 04497 «1.413556 141091 1. 40828 000530 744s 74582 OOO9TH 
cis-1,2- Dimethyleyclopentane CHa 99. 532 04603) «1. 42217 1.41963 1.41709 000508 T7M2 Tht OOO8US 
trans-1,2-Dimethyleylopentane CoHyy V1. SH 4521 141200 «1. 40041 1. 40682 OOO5TS 7424 OOOO 
s-1,3- Dimethyleyclopentane CoHy 91. 725 04525 1.41074 1.40813 1. 40552 oons22 TASS OOOS8Y | 
trans-1,8-Dimethyleyclopentane C Hy 90. 773 (4518 1.40804 «1.40633 «1, 40372 (0522 T4573 000906 i 
Vropyleyclopentane CH 130. 449 O4ASKR 14262 «1 423801 42052 QOO4T4 7HS11 (HON ZZ j 
lsopropycyclopentane CH 128. 419 4915 1.42582 1.423%) 1.42118 000404 THS51 (WOsir I 
1-Methyl-l-ethyleyclopentane CHy 121. 522 (A805 1.42718 1.42476 1, 422384 OO0484 7720) OOS ' 
1- Methyl-cis-2-ethycyclopentane C Hy 128.050 OANYT =1.42983 «1.42805 1. 42457 O04 76 770608 (WINE ' 
1-1-2, Trimethyleyclopentane CH 113. 729 OASIS 1.42208 «1.42051 «1.41804 OOO404 TAYE OOOSRAS b 
1,1,3-Trimethyleyclopentane CsHy 104. 893 (4724 =L4llI9 1.40870) 1. 40621 OO040s TAG5S QOOSHT ' 
*, cis, frans-1,2,4-Trimethyleyclo | 
pentane Cu 116. 731 O4AS2T 1. 41855 1 41612 1. 41369 OOOHSE THA 7a 75MM OOOR4 ' 
trans, cis-1,2,4-Trimethyleyclo 
pert anne CH 1 200 (4738 1 41000 1. 40812 1 4054 COG) 74727 74502 T3873 OOO 4 
' 
ALKYLCYCLOHEXANES ' 
1,1-Dimethyleyclohexane CH 119. 443 0.04920 1.42000 1.42862 1. 42424 0.000476 OL. 78004 O TTATT 0. TT27¢ 0. O0O8 18 ' 
Propyleyclohexane Collys 1M. 724 05200 $1. 48705 1.43478 1.43251 (OED Tuan) TSUTT TRIS (MurTY ; 
lsopropyleyclohexane CoH. 1M. Mi O5210 1.44087) 1.43861) 1. 49635 (WES? sez! TUN Tua MNITTT 
1.3-Trimethyleyelohexane CoH 5. G2 ONY 1.42055 1 25 «1.4249 OOO S60) T7RRS TT4u7 TT) (WITS 
Butyleyclohexane CoH O47 W412 1.44075 1. 43855 1. 4368 TLE TUG1S aN | TOIT wwT42 
lsobut viceyclohexane Coll s21 O5367 1. 43861 1. 43036 1. 43411 OOO450 7W521 TUl4l TAT) TT | 
ec- But vleyclohexane CoH uA O40 1. 4467 144464 1. 44235 QO0438 SI314 SO sO (WMT aU 
But vievelohe vane ( Bt 171. 501 O42 1 44004 1 44475 1. 44242 (wade S1267 SUSuH) SOSLS LL Lees | 
' 
' 
' 
ALKYLBENZENES 
Benzene CH sO 0Ge 0. (4271 1 M2 1. 49702 1. 4047s 0 OO08S OO STUO! OSTSTO OC NONE TLL ' 
Methylbenzene Toluene C-Hs 110. #26 o¢nse 1. 49605 1.40414 1. 49128 (OO58T seid se) SSTTO mus | 
Ethvlbenzene Cul 16. 1st O4SUS 1. 4058s 1 40820 1. 49050 moss SATO? Sate | SAS (WMINTS ; 
2- Dimethyl benzene o X vlene CH 144.411 Ut De Be Le | 15285) (1. NOR OOO] SSL) S75 STIN (WHOIS 4S ' 
Dimethylbenzene n-X ylem CH 130. 104 (M908 «1.49722 1. 4946 1. 49108 OSS SH417 se se HIND? ' 
+- Dimethyvibenzene »X viene CH 138. 351 (M917 1.49582 1. 4932. 1 49054 OOOSRT SOTO See a) so TT wi) ' 
Propyl benzene CH 179. 218 05143 1. 49202 1. 4805) 1. 48702 Tie Lt sud aT S1370 Tiel 
See footnote at end of table. 
' 
' 
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TABLE I Va 


es of 


Pressure 
coefficient 
of boiling 


Boiling 
point at 


Refractive index * np 


Compound Formlua 49 mm point at 760 
He mm He 2° ¢ 
dtdP 
ALKYLBENZENES 
( CimmHg 
lsopropylbenzene (Cumene C,H 152. 392 O5074 1. 49145 
Methyl-2-ethylbenzerne (\H 165. 153 O5163 1. S456 
1- Methyl-3-ethylbenzene CH 161. 305 O5111 1. 49660 
1- Methyl-4-ethylbenzene Coil 161. 989 O5148 1. 49500) 
1.2.3-Trimethylbenzene (Hen 
litene Cyl 176. O84 O5203 1. 51398 
2.4-Trimethylbenzene Pseud 
umene CoHy 169. 351 OSIST 1. SO484 
+ Trimethyvlbenzene Nlesity 
ene Cy 14. 71 OS 100 1. 499387 
i- Butylbenzene Phenyl butane Cyolly ISS. 270) 05358 1. 48979 
I bu vit 7 l 
Phenvl-2-met! propane ( iH 172. 750 OS319 1. 48046 
Butyl nz 2 
Ihe butane Cw 173. 0 531 1. 49020 
Bu t t ne 2 
henyl-2-methylpropane CH 19. 119 05260 1. 49266 
2- Diethylbenzene Coll ISS. 425 05340 1. 50846 
Diethvibenzene ( Ili IAL. 1 O5288 1. 49552 
- Diethyviber ‘ ( His 183. TAZ OASSI 44485 
I i tel hyaer wt he jt tate a itmmospl 


IV. Densities 


The measurements of density were made with 
a density balanee, the assembly, calibration, and 
operation of which has been previously deseribed 


6]. The experimental results on density are 


Individual measurements were 
0.00002 to + 0.00008 2 ml 


given in table 1 

reproducible within 
The accuracy of the tabulated values, including 
the effeet of impurities, is estimated to be + 0.00005 


to + 0.00010 ¢ ml 


V. Boiling Points 


The normal boiling points of these compounds 
were determined as part of another investigation 
on their vapor pressure, and the values reported 
in table 1 for the normal boiling point and the 
pressure coefficient of the boiling point at 1 
atmosphere are taken directly from the report of 
that 
of the 


within 


investigation |7 Individual measurements 
were reproducible 
The ac- 


curacy of the tabulated values of the normal boiling 


normal boring 


0.002 deg to 


point 
0.0038 deg C 


476 





point, including the effeet of impurities, ts 
mated to be 
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Fatigue Characteristics of Electroformed Sheets with 


and without Iron Backing 
By Harry K. Herschman and Carroll Thomas' 


The flexural fatigue properties were determined for electroformed sheets of iron-nicke. 
and iron-nickel-chromium, and of plates composed of electroformed iron nickel sheets bonded 
with solder or a plastic adhesive to open-hearth iron. The fatigue limit of the iron-nickel 
composite was decreased when a thin deposit of chromium was applied to the nickel face. 
Heating the iron-nickel-chromium sheets at 260° C (500° F) before machining improved their 
fatigue properties. Apparent further beneficial effects resulted by heating after machining 
The fatigue limits of the plates bonded with the plastic adhesive were significantly higher than 
those joined by soldering. Two brushed coats of the plastic adhesive applied to the mating 


surfaces of the composite lavers of the plates appeared somewhat more favorable to higher 





fatigue limits than a combination one brushed and one sprayed coat procedure, 


I. Introduction 


Structures or segments of structures frequently 
are prepared by bonding a metal or alloy having 
certain properties to another metal or nonmetallic 
material having wholly different properties. Such 
composites usually will have a combination of de- 
sired characteristics that no single material pos- 
sesses. The joining of the constituent parts may be 
accomplished by various means, such as soldering, 
fusion welding, or adhesive cementing. Many of 
the intaglio printing plates used at the Bureau of 
Engraving and Printing are prepared by bonding 
an eleetroformed face shell composed of a layer of 
iron and a laver of nickel to a backing of open- 
hearth iron. Previously, soldering supplemented 
by spot welding was employed for joining these 
two segments. Because of failures in the solder- 
bonded plates, plastic adhesives have been replac- 
ing this method for securing the required bond be- 
tween the electrolytic shell and the iron back. 

The manner of affixing a plate on the printing- 
press bed causes a slight convexity of the printing 
surface when the plate is not in operation. The 
loading of the plate as it passes under the roll 
While printing tends to flatten the printing surface, 
resulting in a repeated flexing of the plate. It is 
this flexing action that causes failure by fatigue in 





id I man, Electrolytic Section, Bureau of Engraving and Printin 
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some of the constituent layers of the electro- 
lvtically “built-up” shell, which is the portion of 
the plate subjected to the maximum tensile 
stresses. This has occurred particularly in the 
plates bonded with tin-lead solders. A number of 
factors, such as residual stress in the eleetroformed 
shell, the nature of the bonding agent, and brittle 
constituents at the interface of the iron layer of 
the electroplated shell and the solder may influence 
the fatigue properties of these composite plates 
In order to ascertain how these factors affect the 
fatigue characteristies of the electroformed shells 
alone and as bonded to fabricated iron backs, their 
behavior under controlled fatigue stressing Was In- 
vestigated. 

A review of the literature dealing with the 
fatigue properties of electroformed metals and of 
structures bonded with plastic adhesives or solders 
revealed only scant information on these subjects 
A Navy report * on the mechanical properties of 
structures bonded with a plastic adhesive stated 
that a magnesium assembly joined with this ad- 
hesive still carried 100 pereent of the design load 
after 10% 10° eveles, whereas a similar structure 
assembled with rivets failed in) 3% 10° eveles. 
Similar results were reported by Westbrook ° for 


Naval Air Experimental Station, Philadelphia Report (April 1, 1044 
Adhesives for gluing thin metal sheets, F. A. Westbrook, Steel Pro« 
j, 26 (Jan. 194s 
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lap joints secured both by riveting and adhesive 


cementing respectively. Grover and Jackson ‘ 
found that lap joints of magnesium sheet secured 
with “Cycleweld” were stronger both statically 
and in fatigue than joints made by riveting and 
spot welding and were comparable in these res- 
spects with “Heliare’’ welded joints. However, 
data showing the comparative fatigue properties 
of joints bonded with solders and adhesive cements 
respectively, do not uppear to be available. 

Most of the reports discussing the fatigue prop- 
erties of electroplated metals consider them in 
connection with their influence as a coating on 


other metals or alloys. An unpublished report of 


Oberg and Powers?’ indicates that the flexural 
fatigue limit of O.1-in. sheets of electroformed 


43,000 Ib in’. 
F for 2 hr increased the 


nickel is about 
Heating this sheet at 400 


fatigue limit to 46,000 Ib in’. 


“as deposited” 
II. Materials and Specimens 


The test specimens used in the first: phase of 





* The fatigue strength of lap joints it me magnesium sheet alloys, H. J 
er and L. R. Jackson, ASTM Tech. Pub. No. 72 (Symposium on test 
parts and assemblies 
Private communication, T. T. Ot ind J. B. Powers, Wright, Patter 
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this investigation were prepared from electro- 
formed sheets (fig. 1) consisting of a layer each of 
iron and nickel and iron, nickel, and chromium 
Those employed for the second part of the study 
were cut from plates prepared by bonding electro- 
formed shells composed of a laver each of iron 
and nickel to open hearth iron backs (fig. 2). 
The electroformed composite sheets of iron 
nickel were prepared by first depositing a laver of 
nickel about 0.005 thick nickel-faced 
cathode, approximately 13 by 15 in. in area, the 


in. on a 
surface of which was covered with a very thin film 
of flake graphite to facilitate separation of the 
electroformed The then 
transferred to the iron bath, from which an iron 


sheet. cathode was 
layer about 0.1 in. thick was plated on the nicke! 
deposit formed in the first step. The nickel-iron 
composite then was separated from the cathode, 
following which the sheet was ground on the iron 
surface to a finished thickness of approximately 
0.075 in. 
composites, the nickel and tron were plated as 
Then the nickel surface was 


In preparing the iron-nickel-chromium 


previously described, 
electrolytically cleaned with a solution of 45 g¢ of 
sodium hydroxide and 25 g of trisodium phosphate 
per liter of water, scrubbed with pumice moistened 


with 20 percent of hydrochloric acid and finally 








Ficgure | Sections of fatique specimens of electroformed sheet shou ing thickness of the nickel and chromium 
layers respectively. 
The lave f copper and nickel indicated at the top of the micrographs were applied to the metallographic specimens to prevent rounding of the adia 
1 {the specimens during p hit Number designations of micrographs correspond to sheet numbers in table 1 Etched in 5 percent picra Ma 
ation Th 
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the sheet in the chromium plating bath. (3) Chromium bath 
The conditions of plating and the compositions CrO 250 g liter 
for the different electroplating baths employed H.SO, 2.5 g liter 
were as follows: Current density =27 amp square decimeter; : 
(1) Nickel bath Temperature = 50° C; 
NiSO,.7H.O ;, 200 ¢ liter Time of plating =40 min. 
NiClL.6HLO aoee 10 ¢ liter The thickness of each component laver of the 
HBO 20 ¢ liter clectroformed composites, tested as such, are listed ' 





Current density = 2 amp square decimeter; in table 1. The test specimens prepared from this 
pH=5.1; sheet material were of the form and dimensions 


: 
Fiat RE 2 Sections of fatique specimens of printing plate com posed ofan electroformed portion ofa layer each of iron and 
nickel bonded to a backing of open hearth tron. 
Number designations of micrographs correspond to plate numbers in table 2 Magnification «4 
rinsed with water immediately prior to immersing Time of plating=36 hr. 


: 

Temperature =40° C; shown in figure 3. A number of the specimens 

Time of plating =6 hr. prepared from the composite sheets 2) and 3 i 

2) Iron bath (table 1) (Numbers 2B and 3A, table 1) were 

FeCl, 4H.0 450 & liter tinned at 260° C (500° F) on the iron face in some | 

CaCl, 200 ¢ liter cases and others (number 2A, table 1) were 

Current density =7 amp square decimeter; heated at this temperature without tinning. The ; 
pHi=1.0; edges of the specimens within their gage length 
Temperature =95° C; were polished with 3,0 fine emery paper to remove 
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surface irregularities that might serve as centers of 
stress concentration and thereby act as nuclei for 


fatigue failures 


Tarte 1 Fatique mits of electroformed composite sheets 
on-nickel and tron-nickel-chromium 
Thick 
Fatigue 
t Elec lin 
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Parte 2 Results of flerure fatique tests on bonded plates 
] | 


The eleetroformed shells for the bonded plates 
were prepared in a manner essentially similar to 
that used for the electrolytic sheets deseribed 
above In every case they consisted of a laver 
each of iron and Tile kel ol the respective thiek- 
nesses indicated in figure 2. These composites 
were ground on the iron face with a Blanchard 
grinder, and this surface then was cleaned by 
Finally 


brushing with lime, pumice, and water 


it was washed with a 10-percent aqueous solution 
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Figure 3. Design of specimen used for determining 


fexcural fat que prope ties of electroformed sheets 


of hydrochloric acid. The bonding surfaces 
the open hearth iron backs were similarly ground 
cleaned, and treated with acid 

Soldering the electrolytic shell to the iron ba 
was accomplished by fluxing the contacting iro: 
surfaces with zine chloride, sandwiching a 0.005 
in.-thick sheet of 50-50 tin-lead solder betwee 
these surfaces, pressing them together at a ten 
perature of 260° C (500° F 
80 Ib/in- 


heated plattens and cooling under pressure of 8 


under a pressure 


in a press equipped with eleetrically 


Ib/in.? in a water-cooled press 

When the joints were made by adhesive cemet 
ing, the mating surfaces first were washed wit! 
an abrasive soap and water and then clean: 


After 


a coat of plastic adhesive was applir 


with steel wool thoroughly drying thy 
surfaces, 
with a brush to each of them. Following a dryvir 
period of 30 min, a second coat of adhesive w: 
applied, in some cases by brushing and in othe 
by spraving, which also was allowed to dry in 
for 30 min. These coatings then were dried 
85° C (185° F) for 25 min, after which the coate: 
surfaces were mated and pressed together tine 
a stress of approximately 80 Ib in’ at a tempel 
ture of 160° C (325° F) for 30 min 

The plates bonded by these methods had 
nickel face and an iron back. Surface metal w 
removed from the iron backs by grinding with 
Blanchard surface grinder until the plates we 
thick 


form and dimensions shown in figure 4 were « 


approximately QO.2 in Specimens of t 


from these plates. In all cases, the edges of t! 


specimens within the gage length were polish 


with fine (3.0 emery) abrasive paper to remo 


surface irregularities in the same manner and 
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FIGURE } Design oj} specimen used for determining the 
flexural fatigue properties of composite plates of electro- 


formed nickel-tron shells bonded to open hearth rron backs 


the same purposes as described for the electro- 


formed sheet specimens. 


III. Method of Test 


The flexural fatigue tests were made on Krouse 
plate fatigue machines operating at 1,750 rpm 
fig. 5). The electroformed sheet and composite 
plate specimens in every case were adjusted in 
the machine so that the range of stress applied to 
the specimen face (nickel or chromium surface 
see table 1) was from zero to maximum tension 
and the back (iron surface) was in zero to maxi- 





Fieure 5 Vachine used for making flexural fatique test 
he nag manne i? hicl Specimen ts suppo ted 
pecime « for holding specimen in rigid position; w=secrew for 
ijusting height of vise v; ¢=crank-arm for actuating specimen i=connect 
¢ pin for attaching arm ¢ to specimen s; ¢ uijustable eccentric for varying 
ow of connecting arm ¢; r=reflecting mirror used for ascertaining spe 
flection in terms of stress as shown in figure ¢ 
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mum compression. This stress cycle was used 
to simulate that to which a printing plate is sub- 
jected on the press. 

The specimens were deflected to produce the 
maximum flexural stress by adjusting the eccentric 
except that the 


(fig. 5, e) in the usual manner ® 


method used for measuring the deflection (fig hy 
was somewhat different from that generally em- 
ploved with the Krouse flexure machine. Essen- 
tially it consisted of a flat mirror (7, fig. 5), attached 
to the adapter securing the specimen to the crank 
arm, which reflected a horizontally collimated 
beam of light to a 12-ft-long seale placed about 
20 ft from the mirror (fig. 6). This permitted the 
reproducibility of any desired maximum stress to 
within 1 percent. 

Each testing machine was equipped with a 
microswiteh, which was automatically operated to 
open the motor circuit and thus stop the machine 
when a test specimen ruptured completely. A 
counter attached to the machine recorded the 
number of flexural stress eveles to which a specimen 
was subjected during test In the case of some of 
the plate specimens, particularly those bonded by 
soldering, the electrolytic shell and the bond 
frequently failed first, whereas the back remained 
Intact Under this condition the stress in the iron 


back was decreased below its faticue limit, and 


y, 
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the machine would run until stopped manually. 
Therefore, unless an operator Was in constant 
attendance, a record of the number of cycles for 
failure of the shell was not obtained. In order to 
provide for this contingency, a special circuit 
comprising a relay in series with a number 38 B 
and Svage silk covered enamel copper wire secured 
with Duco cement to the upper surface of the 
specimen (fig. 7) was installed to interrupt the 
test when rupture of the shell occurred before 
complete failure of the back 
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r < 
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Figure 7 Schematic diagram of the relay sustem used to 
liscontinue the flerural fatique test hen the electroformed 


nickel-iron lane A hite ection failed while the fabri- 


cated iron backing B (cross-hatched) remained intact 

When the electroformed lave A fractured, the copper wire W (grounded 
from G through the mact ruptured, opening the relay R and stopping 
t) } 


In the present study the fatigue test was em- 
ploved primarily as a means of obtaining an ap- 
proximate relative measure of damage by flexing 
associated with the plating, thermal and bonding 
Because of this and the 
limited number of specimens available for test, 


treatments emploved 


the values selected for the fatigue limit were based 
on a stress, 500 to 1,000 [bin below the lowest 
stress resulting in failure, provided at least two 
runs of not less than 10 million eveles had been 


made at this stress 


IV. Results and Discussion 


The S-N (stress-cvele) curves for each of the 
electro-formed materials tested are shown in 
figures S to 12, inelusive. The fatigue limits 
indicated in these figures are listed in table 1 
It is apparent from a comparison of the curves in 
figures S and 9 that the fatigue properties of an 
iron-nickel composite sheet were affected adversely 
by the application of a thin (0.0003 in.) electro- 
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STRESS CYCLES 


Fiaure 8&8 Stress-cycle (S-N) graph based on results 
flexural fatique lesis made on speciemens of sheet No 


table 1). 


deposited coating of chromium (fig. 1, b lHlows 
ever, this influence of the chromium layer wa- 
somewhat beneficially modified by thermal treat 
ment (heated at 260° C (500° F) for 5 min 

(figs. 9and 10). Maximum improvement appeared 
to be obtained when the chromium plated tron 
nickel sheet was heated (260° C for 5 min) afte: 




















machining. In the cases of the solder-coated 
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STRESS CYCLES 
Ficure 9 Ntress-cycle S-N) graph based on results 


fleaural fatique tests made on specimens of sheets No 
2A, and 2B (table 1 


us electroplated (sheet 2); X, tinned 5 min at 500° F after mac 
ing (sheet 2A); @, heated 5 min at 500° F after machining (sheet 2B 
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specimens, any deleterious effects that may have 
been imparted by the iron-tin compound formed 
at the iron-solder interface, apparently were over- 
shadowed by the improvement resulting from the 
thermal effects (figs. 9 and 10 

The results of the tests on the specimens 
machined from the electroformed sheet 4 of iron- 
nickel-chromium after solder coating (fig. 1, d) at 
260° C (500° F) (at temperature 5 min) showed 
considerably lower fatigue properties (fig. 11) than 
aun of the eleetroformed sheets tested. It is 
noteworthy that the fatigue properties were 
significantly higher for the specimens cut from 
sheets 2B and 3A (table 1, also fig. 1, b, and ¢ 
having a 0.0003-in. thickness of chromium plate 
than for those prepared from sheet 4 (table 1 and 
fig. 1, d) plated with a 0.0006-in. thickness of 
chromium. Recent work done at the National 
Bureau of Standards * showed in certain cases that 
a plate thickness of 0.001 in. of chromium reduced 
the endurance properties of steel more than an 
appreciably thinner electrodeposit of this metal. 
This suggests that the lower fatigue limits (table 1) 
obtained on specimens of sheet 4 as compared 
with those machined from sheet 3 were attribut- 
able, at least in part, to the thicker deposit of 
chromium on sheet 4 























ffect of chromium plating on the endurance limit of steels used 
H. L. Logan, J. Research NBS 43, 101 (1949) RP2O11 
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STRESS CYCLES 
Ficure 10.—Stress-cycle (S-N) graph based on results of 
flexural fatique tests made on specimens of sheets No. 3 
and 3A (table 1 


Heated for 5 min at 500° F before machining (sheet 3); @, heated for 
it 0° F, machined and then tinned for 5 n.in at 400° I sheet 3A 
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STRESS CYCLES 
Fieure 11 Stress-cycle (S—N graph hased on the results 
of fleaural fatique tests made on specimens of sheet No 
table 1 


The fatigue limit of sheet 5 (fig. 1, e) formed by 
continuous electrodepositing (no detreeing) was 


somewhat lower than for material electroformed 
» 


similarly (sheet 2, table 1, fig. 1, b) exeept for 


intermittent removal from the plating bath for 
detreeing However, the adverse effect on the 
fatigue limit by thicker deposits of chromium 
previously discussed is considered to have been a 
greater influencing factor than the omission of the 
detreeing treatment, 
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STRESS CYCLES 
Figure 12 Stress-cycle (SoN graph hased on the results 


of flerural fatique tests made on 8 PECLMENS of sheet Ne 


table 7 
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STRESS CYCLES 
Fiaure 13 Stress-cycle (S-N) graph based on the results oJ 
fatique test made on pecimens oft printing plate 
composed of an electroformed shell of a layer each of iron 
bonded to a backing of open hearth iron 


and nickel (hq 


th two brushed coats of plastic adhesive 


The results of the tests on specimens prepared 


from plates composed of electroformed sheets of 


iron-nickel bonded to open hearth iron backs are 
shown in figures 13 to 15, inelusive. It will be 


noted that the solder-bonded plate (fig 15) had 
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STRESS CYCLES 


Fiaure 14, Stre -cucle SN) @ apl hased on the result 


of flexural fatique test made on specimens of printing 
ate composed of an eelctroformed shell of a layer each of 
on and nickel fig. handed to a hael ing of open hearti 


on ith wt tiah shed coat and fhen a sprayed coat oft 


lastic adhesive 
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an appreciably lower fatigue limit than either « 
the composites bonded with the plastic adhesiy 
13 and 14). 
that were bonded with two brushed coats of th 


(figs. Furthermore, the specimet 
adhesive showed superior fatigue properties t 
those joined after applying first a brushed coa 
and then one by spraying. In nearly every cas 
the failure of the soldered specimens was accon 
panied by a failure of the bond (fig. 16, a). Th 
plastic bonded specimens rarely showed any in 
dication of this type of failure, complete ruptur 
occurring without separation of the componen 
electroformed and fabricated iron lavers as show: 
in figure 16, b. 

A factor that in the greatest measure may by 
responsible for the adherence of the layers bonde: 
thos 
soldered is the degree of continuity of the bond 


with the plastic and the separation of 
Thus, in every case examined the plastie bond 
was continuous (fig. 17, b) and in the case of 
separated soldered joints, there was evidence o 
discontinuities or porosity (fig. 17, a) 
laver. Another factor that 
this result is the fact that the modulus of elasticity 


may contribute t 


of the iron-tin compound at the interface of th 
solder and iron is higher than that of the plastiv 
adhesive. 

During the flexing action of the specimens, the: 
Was an appreciable stress applied to the bond 











laver. If the bond layer is of a porous nature, such 
as illustrated in figure 17, a, then its fatigu 
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strength might be sufficiently low to permit rup- 
ture during flexing with consequent separation of 
the two constituent lavers. Despite this separation 
the lavers remained firmly cemented togethe ut 
their ends. This permitted a slight elongation and 
bowing of the eleetroformed laver fig. 18), thus 
Hnposing a Compressive stress on its bottom surface 
and tensile stress on its upper surface when the 
specimen was in the initially set zero stress posi- 
tion \s the specimen passed through i complete 
flexing evele, the bottom face of the cles troformed 
laver being above the center plane of the specimen 
parallel to its broad face (fig. 18) was subjected to 
a reversal of stress from compression to tension to 
compression This suggests that since the maxi- 
mum tensile stress imposed on the bottom face of 
the electroformed layer was approximately thi 
sume before and after the solder layer failed, the 
stress range for this bottom face increased after 
bond failure This factor, coupled with stress 
raisers on the bottom face associated with the 
voids in the solder lave fig. 17, a) and the corro- 
sion effects of the excess soldering flux residing in 
the voids, may have caused fatigue failure of the 
clectroformed lave 

Although the bonded printing plates commonly 


are chromium-plated on the nickel face, none o 
the bonded plates tested were thus plated The 
results obtained on the electrolytic shells (figs. 8 
to 12, inclusive) clearly indicate that the endurance: 
properties of chromium-plated bonded — plates 
would be somewhat lower than those shown for 
the bonded specimens tested (figs. 13, 14, and 15 


However, on the basis of available service data, it 
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is not expected that the chromium plate would 
have altered the relative results as influenced by 


the type and method of bond employed. 
V. Summary 


\ laboratory study Was made to ascertain the 
flexural fatigue properties of electroformed = iron- 
nickel and iron-nickel-chromium composite sheets 


and of plates composed of electroformed iron- 


3. The application by brushing of the two coats 
of plastic adhesive to the mating surfaces of the 
composite lavers of the plate resulted in a higher 
fatigue limit than the one brush and one spray 
coat procedure, 

$. In nearly all cases, the solder bond failed 
either just prior to or simultaneously with the 
failure of the electroformed layer; the fabricated 
iron back remained intact under these cireum- 


stances. The number of plastic bond failures was 
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Fiatre IS schematl fiagqran ndicating the 


ormed shell due to elongation of the she 


the electrotormed she 0 


nickel sheets bonded to open hearth iron The 
electroformed sheets in every cause consisted of one 
The bond- 


ing agents used in preparing the plates were 50-50 


laver of each of the constituent metals 


tin-lead solder, and a plastic adhesive applied ia 
some cases as two brushed coats and in others as 
one brushed and one sprayed cout respectively to 
each of the two mating surfaces. On the basis of 
the results of this investigation, it is coneluded 
that 

|. The fatigue limit of an iron-nickel composite 
is adversely affected by the application of a thin 
electrodeposited coating of chromium Heating 
260° C (500° F 


such tri-laver composites at 


caused some increase of its fatigue limit. The most 
pronounced improvement resulted from heating at 
this temperature after machining. 

2. The fatigue limits of the plates bonded with 


the plastic adhesive were significantly higher than 
those of the plates joined by soldering 








howind eflec 
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{ produced in the electro- 


probab 
ing flerina atte fa ire of the bond but before 


on bach 


negligible: inh nearly all CuSCS the Iron lave re 
mained secured to the electroformed laver and 
they failed as a unit 

5. The lower fatigue strength of the solder 
bonded specimens Was attributed to the britth 
Iron-tin laver at the interface of the solder and 
iron, porosity or discontinuity of the solder lave: 
and corrosion caused by residual flux residing in 
the voids of the solder 

6. The results of this study are consistent witl 
the service data of printing plates compiled at the 
Bureau of Engraving and Printing over a period 
of vears in that the average life of the plastiv 
bonded plates is considerably greater than the 
uverage life of plates whose component lavers were 


joined by solder nye. 


Wasuincron, May 27, 1949 
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Interfacial Instability and Mixing in Stratified Flows 
By Garbis H. Keulegan 


When a liquid flows over a body of heavier liquid, with which it is miscible, an interface 


may be formed at which there is a sharp discontinuity of densities but not necessarily of 


velocities As the relative velocity is increased, waves are formed at the interface, and at 


a critical velocity, mixing begins by the formation of eddies periodically ejected from the 


crests of the waves The numerical values of a criterion for mixing, whose form was derived 


from theoretical considerations, were determined experimentally Experimental data were 


also obtained on the relation of the amount of mixing to the relative velocity of the liquids 


at velocities higher than the critical 


I. Introduction 


The intrusion of salt water in rivers that com- 
munieate with the seas and the silting of reservoirs 
receiving muddy waters are the two most impor- 
tant examples of stratified flow where the question 
of interfacial behavior assumes a basic signifi- 
cance. Depending on the conditions of approach, 
two possible forms of the interface may be noted 
In one form, the interface may be identified as the 
dividing surface of two layers of liquid with dif- 
ferent densities, the surface being one of sharp 
discontinuity of densities but not necessarily of 
velocities. Ordinarily the interface of this type 
is the locale of internal waves if the difference in 
velocities at points on Opposite sides of the inter- 
face and at some distance from it is large. If 
mixing is present, it is in the form of eddies that 
are pet odically ejected from the crests of the 
waves into the current that has the greater ve- 
locity. In the other form, the interface ts a laver 
of transition between two currents. Both the den- 
sities and the velocities change uniformly in the 
laver that has measurable thiekness. If any 
mixing is present, it is associated with the mo- 
mentum exchange of turbulence, and the regular 
pattern of internal waves is absent. Thus a 
study of interfaces may be confined to one or the 
other of these forms. 

In the present investigation, which is restricted 
to an interface with sharp discontinuity of densi- 
ties, one may visualize the following distinct 
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problems. First, what quantities are to be con- 
sidered in determining a criterion for the incidence 


) 


of mixing? Second, with current velocities exceed- 


ing the critical velocity of the criterion, how is 
the amount of mixing related to the velocity? 
Third, what is the exact mechanism of flow that 
produces the mixing? 

About a decade or so ago at the National 
Hydraulic Laboratory, experiments were carried 
out to study each of these three questions. After 
having studied the first and the second of them 
with some degree of satisfaction with the results 
obtained, we were about ready to enter into the 
study of the third, which obviously is beset with 
considerable experimental difficulties and theoret- 
ieal uncertainties. Unfortunately at that) time 
it became necessary to suspend the hydraulic 
work in the laboratory, and since then we have 
had no opportunity to return to our original plans 
Meanwhile considerable interest in these matters 
has developed. Since numerous inquiries have 
been received about the progress already made, 
it was decided to make these results known not- 


withstanding their incompleteness. 


Il. Determination of the Critical Velocity 
of Mixing 

Observations on the instability. of interfaces 

were made for the condition where the upper, less 

dense liquid was flowing over a heavier liquid at 

rest inva pool, The essentials of the apparatus 
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and the arrangements are shown in figure 1. The 
apparatus consisted of a closed flume with vertical 
vlass walls and was limited at the entrance by a 
constant level tank and at the exit by a discharge 
tank The large reservoir of the circulating 
liquid, which was the lighter liquid, was connected 


to the constant level tank through a pump. To 





The flowing current caused the lower liquid t 


be dragged out in the downstream direction lt 
the small and the medium flumes, the escape of 
the lower liquid was prevented by tilting the 
flumes by means of screw-jack supports. This 
method was not practicable for the large thumy 


Instead, the flume was fixed in a moderately 
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prevent the transmission of pump vibrations to 
the flume frame, the connection to the constant 
level tank was made through a flexible rubber 
tubing. The return flow passed from the cis 
charge tunk into a metering tank of the orifice 
type, which in turn discharged its contents con 
tinually into the reservoir of light liquid below 
The volume of the circulating liquid during a test 
remained the same, and since the reservoir volume 
Was many tithes larger than the volume of the 
flume, mixing in the flume during a run had very 
little effeet on the density of the upper liquid 
The loss of the heavier liqguid from the pool during 
tests Was made up at intervals by drawing from a 
reservoll placed uta higher clevation 

For the purpose of examining the effect” of 
Revnolds number, three flumes of various sizes 
were employed. They were rectangular in cross 
section, With the greater dimension vertical, with 
cross sectional dimensions of 2 by fem, 4 by S 
em, and 11.3 by 28.5 em. The length in each 
ease Was about 25 times the total depth The 
flumes were operated full, keeping the depths of 
the two lavers equal or nearly equal. Accord 
ingly, the hydraulic radii associated with the cur- 
rents were 0.50, 1.00, and 3.27 em The fumes 
will be referred to as the small, the medium, and 


the large flume 
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When the velocities of the 


current were too high, the lower liquid that was 


inclined position 


dragged out Was prevented from enterines the 
lower reservoir by drawing it off through a trap 
at the exit end. The idea involved in the metho 
is shown in figure 2 

The determination of the critical velocity; tha 
is, the velocity at which the mixing of the liquid 
of the lower laver with the upper commenced, o 


was about to commence, was done visually by 

















noting the condition of the interface At low 
velocities a smooth interface was discernible ever 
When the lower liquid was not colored. At th 
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initial stage, observations showed that at the 
interface, and in a finite band on the two sides of 
the interface, the flow was laminar even when the 
central regions of the upper current were eddying 
or turbulent. When the velocity of the current 
was increased to some definite value depending 
on the densities and the viscosities of the two 
liquids, the interface appeared to be covered with 
ridges extending from one glass wall to the other. 
These ridges moved progressively downstream, 
and their crest lines were parabolic in shape, the 
crest heights being greatest at the center. With 
a slight increase in the value of the current velec- 
ity, waves appeared in a new configuration. The 
interface now was covered with waves of shorter 
crest lines. Inthe tests with the large flume these 
crest lines were one-half or one-third of the channel 
width. The appearance of the waves and their 
sequence was regular. The movement was pro- 
gressive. The predominant wavelength remained 
constant, or nearly so, throughout the channel 
length. With further increase of the velocity of 
the current, the waves became sharp-crested, and 
mixing commenced, At the beginning the mixing 
was slight, and it augmented with subsequent 
increases in velocity. [It appeared that increased 
velocities did not substantially affect the wave- 
lengths of the interfacial waves, although there 
was a tendency for the lengths to decrease slightly. 
The most marked influence of the increasing veloc- 
ities Was to augment the height of the sharp- 
crested waves, the size of the eddies emanating 
from the crests, and the frequency of departure of 
the eddies. Practically, the mechanism of mixing 
Was an ejection of eddies from the crests into the 
moving current, as shown diagrammatically in 
figure 3. 

Now, what constitutes the eritical velocity of 
mixing isa matter difficult to decide. The method 
that we employed is as follows: In the record book 
was entered a deseription of the appearance of the 
interface, together with a qualitative statement 
regarding the degree of mixing for a measured 
velocity of the current. The interface, for exam- 
ple, was described as being smooth, or having para- 
bolic waves, or mixed waves, or sharp crests, or 
agitated crests. The mixing was described as 
absent, or slight, or moderate, or eXCeSSIVe., We 
note the velocity Uy, the maximum velocity record- 
ed in the book as corresponding to a smooth inter- 
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face and to no mixing. Next we take the velocity 
U\, which is the average of all the velocities asso- 
ciated with slight mixing. The mean value of 0) 
and (7), is chosen to represent the critical velocity 
of mixing. It is denoted by U,. The average 
value of the velocities associated with moderate 
mixing are denoted by Uy, and the average value 
above which there was always excessive mixing 
as (%,. 

The test liquids were tap water and sugar solu- 
tions. Four series of tests were made. In three 
of the series the upper liquid consisted of tap water, 
and tests were made in all three flumes. The pur- 
pose of these tests was to determine the effect of 
the Reynolds number of the upper current on the 
incidence of mixing. The fourth series of tests 
was intended to determine the effect of \ ISCOSILV 
on mixing. These tests were made with the medi- 
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Figure 3. Diagram of process of miring 


um flume. Both the upper liquid and the lower 
liquid were sugar solutions. The fractional differ- 
ence of the densities was kept practically constant, 
but the viscosity of the upper liquid was varied 
The velocity data from the tests of the four series 
are given in tables 1, 2,3, and 4. 


TABLE 1 Miring tests in the small flume 
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Tair 2. Mixing tests in the medium flum III. A Theoretical Basis for the Dimension 
Hydraulic radius= 1.00 en less Parameter of the Critical Velocity 


If it were not for the inertia and the restraini: 

effect of the upper current on the motion of th 

ms m cand arr hae gd baer interfacial internal waves, we would introdu 

D-2 seo | 1.000 | .eune | cree | aa ; r ; . directly the criterion for wind-generated wav: 
D OCS | 505s | -C0es | .Cre | ae wae = : derived by Jeffreys as the desired criterion 1.2 
Jeffreys has shown that at the moment « 


9 ‘ 6.4 2 developing wind waves 


D9 ONS »” 2 ie , ¢ > 54 O04 & 80 
es of Frapersal Rapnid figrenel Mash Baie sean Bieacnad lbghns where (7 is the velocity of wind over water that 
D-12 6 LISS 0276 ONO 8 “9 Lm 9 still initially, » is the kinematic viscosity o 
D » te a . : oo ee ~~ Water, Ap Is the difference between the densitix 
Dou 20 02 1046 978 Hos 12. of water and air, p; is the density of air, S is : 
< a aa pan : “ ~ , op wen , ~ numerical constant, and g is the acceleration o 
gravity. 

Jeffrey’s method for the analysis of the win: | 
TABLE 3 Miring tests in the medium flume velocity criterion shows a marked departur 
Hydrau 00 en from the well-known derivation due to Kelvir 
The hypothesis that the flow in air is irrotationa 
is abandoned, and instead it is supposed that th 
motion of air ts such as to introduce a pressu 


18 OOo 1. 268 0.008 su 2 912 1 Variation over the water surface of the magnitud: 
wm) 02 SS HOE MY KS p—Sp,l"dh/de, | 


where S is a numerical constant referred to as th |, 
‘ ; sheltering coefhicient, [ is the velocity of the a 1 
a ' ; '. 6.8 3| s00!) 980! 11.38 over the crest, A is the surface elevation of thi is 
wave above the undisturbed surface, and wv ts i r 
the direction of wind motion ti 
Poeee 4 Aficine teste ja the fora: Gum Jeffreys has concluded from observations o re 
the initial formation of wind waves ina river and 
the waters of a large pond that the velocity of 7 
wind just strong enough to raise waves ts 11 0 
em sec, the average of threeobservations. Thecor | 
responding value of S is 0.274. Substituting thes di 
. values in eq 1, we get for the criterion, writing ti 
) ( vo *) / : t. 
' ' , m a oR , fA p 
the value 


| 7.98 @—0.215 


The problem of interfacial waves between tw a 


so) 629) #90 moving lavers is a question of small oscillation my 
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about a steady state of motion. Let the velocity 
of the upper laver be 0’, and of the lower layer zero, 
Let the densities be  p, and 
po(—Ap+p,). The wavelength and the period 
being Aand 7, the internal waves may be assumed 


corresponding 


to be given by 
h—a cos (kr—ot). (5) 


where K=22 \ and o=22/7. Supposing that the 
motions are irrotational everywhere, that the 
wavelengths are small in comparison with the 
depths of the layers, and that the wave height is 
small with respect to the wavelength, analysis 


vives for the velocity w of the waves 


o U/ l q Ap _ l v9 3 . 
Ww I: 2 a E k p1 4 l ] ’ (t)) 


provided that the square of App; is negligible with 
respect to unity [3]. In stable waves ¢ is real, and 
thus 
/ 29 Ap 
—k »p 


As was noted previously, the internal interfacial 
waves just prior to the setting-in of mixing are 
regular in form and move progressively down- 
stream without changing their form. The regu- 
laritv and the permaneney of form indicate that 
there isa mechanism by which the upper current 
is furnishing energy to the wave motion of the 
lower liquid just sufficient to counter the dissipa- 
tion due to effects. three 
conceivable methods of furnishing the energy. 


7 

viscous There are 

The energy supply may be effected through the 
1 

normal stresses p.., or the tangential stresses p,., 

Since for 

the ultimate purpose of the analysis, which is the 


or by the combination of these stresses. 


derivation of dimensionless parameters, the selee- 
tion is not important, it will be supposed that the 
transfer is made by the normal stresses p.., and the 
vanish. 


The dependence of pP 


tangential stresses Pre 
on the wave form must 
be known. The dissipation of waves through 
viscous forces has a direct bearing on the relation, 
and it is necessary that the ordinary equations of 
motion involving the viscosity terms be considered. 
These equations, neglecting the second-order terms, 
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It is implied that 2 is drawn vertically upwards. 
It will also be supposed that the origin is in the 
undisturbed level of the interface, and that 7 is 
The 


quantities « and w denote the particle velocity 


taken in the direction of wave motion. 


components along the axes s and 2, respectively. 
It has been shown by Basset [3], that the above 
set of equations is satisfied by 


Od Oy _ 0 , oy 
7 Or Of ’ “O27 Or’ (tt) 
and 
) O> _ ge (12) 
/ p or Y-, _ 
provided that 
; wv, F 
V7o=0, ay TVA. (13) 


Basset [3] and later Lamb [4] have applied these 
to the determination of the damping of the surface 
waves of a homogeneous layer. In reworking the 
solutions, we wish to find the condition that corre- 
sponds to no damping. It is to be noted that 
only in the event that y becomes negligible may 
the flow be regarded as irrotational. 

It seems that the above equations and the solu- 
tion derived from them can be applied directly to 
the case of internal waves after resorting to the 
following artifice. The presence of the upper cur- 
rent is imagined to be evidenced in two Ways only. 
First, the dynamical action of the current consists 
of the normal stresses p,,. Second, the intensity 
of gravity in the lower liquid is reduced from g to 
q’, so that 


q’ hw. (1 4 
Accordingly, in the place of eq 12 we will write 


F Oo ? 1h 
/p ome ¢ Pa ( » 
- or . 
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For the case of infinite depth and liquid at rest 


corresponding to the wave profile, 
h—a cos (kr—ot), 5) 


correspond the solutions of eq 11 


= Ae“ - 16 
and 

y= Cee 4" la 
with 

m k- av. IS 


If the case of cellular waves be excluded, m de- 
notes that root which has its real part positive 


oan 
Thus, 


u the Ne me é — 19 


and 
ihCe™) ett +o. 4) 


In this solution we are supposing that there is 
no bodily movement of the lower laver as a whole, 


so that the surface kinematical condition ts 
Oh /Ot = w, (), 21 
and this leads to 


h=— Ail \e ma 22 


By Stokes’ rule, the relations between the 


stresses and the dilatations are 


P p ak Qu ’ ~*) 
and 
ow . Ou ” 
p u(5-+5- } 24 


As regards the normal stresses, it is appropriate 
to suppose, since a regular progression of waves is 


imagined, that 


Pe BAe** tes, 25 


It being the convention io regard pressures as 
negative tensions, 8 is a positive quantity. The 
dynamical surface conditions, from eq 15, 19, 20 
22, 25, 24, and 25 are 


IvkmayC=—0, (26 


a vk a yk Da | gk 


and 
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The elimination of the quantities land © leads ¢ 


a trk?a gk a8 — 2vk?B4+- 47k ty-k'm —O 


DK 


where v is small, the terms containing its squar 


mav be neglected, and 


a ty a gk ap Dvh re ) 29 
This may be written 
a $(vk?—B/4)a qk vk? 3 


The condition that the waves are neither dampx 


nor growing in height requires that 


Thus, eq 30 reduces to 


a ak, 2 
or 
a or 
and substituting the latter in eq 22, since of 
l 
h A—iC\« 


Kquation 27 shows that the ratio ().1 is a sma 


quantity when vis small. Neglecting ©, 


Putting 


h a cos (ka at 
To this corresponds the surface pressure variat! 
p tp.vk?wa sin (h at 
and therefore the pressure is smaller on the f 
ward face of nwave and vreater on the baa kwy 


face. as is to be expected 
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Here, following Jeffreys, we introduce the 


supposition that 


p Sp, (l —w)*dh/dy, 34) 
that is, 
p Sp U'—w)? ka sin (kx ot). STH) 


The implication is that for the maintenance and 
for the growth of waves U7 is larger than w; that 
is, the current velocity is larger than w. When U’ 
is smaller than w, the waves are damped. The 
latter statement is due to Kelvin [5]. 

Comparing eq 33 and 35 with each other and 
ignoring the small term involving Ap/p,. and 


writing v, for: 


This is the relation to be satisfied between the 
quantities w, (7, and & at the instant that the 
internal waves are about to increase in height 
Returning to the inequality, eq 7, it may be sup- 
posed that the wavelengths have reached the 


greatest values consistent with stability; that is 


2q Ap 
” p - 


Eliminating & between the latter two equations, 


s voJ Ap, pi, (3S 


and sinee under these conditions 


see eq 6) we finally obtain 


16 
S vogAp p 10) 


This may be written as 
0 vogAp p { > +] 
and is the desired form for the criterion of mixing. 


IV. Experimental Value of the Criterion of 
Mixing 

To generalize the result of the theoretical anal- 

ysis made above and in particular to examine the 

effect of the Reynolds number of the current on 
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mixing, the data of the tests will be considered in 
the parametric forms 


O= (rg Ap p;) ‘2 
and 
0—0(U Rv 13 


In these, p; is the density of the upper flowing 
liquid, »; and v. are the kinematic viscosities of the 
upper and the lower liquids, respectively, Ap is 
the difference in the densities of the two liquids, 
(’. is the critical velocity of mixing, and 7? ts the 
hydraulic radius of the cross section of the flowing 
current. 

However, prior to the consideration of the mag- 
nitude of the criterion it will be helpful to present 
the evidence of the effect of kinematic VISCOSITY 
on the criterion itself. The criterion contains the 
quantities vy, and Ap p, each raised to its one-third 
power. The effect of the viscosity will be seen 
best and also most convincingly if a series of tests 
are conducted keeping Ap p, constant or nearly 
constant and varving the kinematic viscosity », 
over a large range. Most of the runs of the J) 
series In the medium flume fulfill this condition. 
The test data are given in table 2.) The plotting 
shown in figure 4 is from this table. The critical 





Uc - cm per sec 











velocity of mixing is treated as a function of the 
kinematic viscosity. of the lower liquid. The 
great majority. of the points aline themselves 
along a straight line having an inclination of one 
in three. Since the plotting is done logarithmi- 
cally, this at least confirms the requirement of the 
theory. 

As regards the criterion and the effect of Revy- 
nolds number on it, the values of the parameter 
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© from all the tests are shown in figure 5. One 
vlance at the figure suggests that the criterion 


has one constant value below a Reynolds number 


of 450 and another constant value above it: 


6=0.127, R,< 450; (44) 
ana 
6—0.178, FR. 150 $5) 


Significantly, 7, 
terion that separates the regime of turbulent flow 


$50 is also the Reynolds cri- 


from the regime of laminar flow 

Considering the distribution of the © points 
more carefully for Reynolds numbers above 450, 
some systematic variations from the mean value, 
O.178, are detectable. This would mean that 
there are other factors having a secondary influence 


on the value of the mixing criterion. Such factors 
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may very well be (a) the surface tension between 
the two liquids, (b) the difference between the 
kinematic viscosities of the two liquids, and (¢ 
the relative densities of the two liquids. It 
would be a futile effort to speculate on the effect 
of these factors separately. Since in these present 
tests the solutions were prepared from the same 
solvent and solute, one may say that all these 
Thus 
it would be appropriate and also sufficient to 
For that 
purpose the graph in figure 6 has been prepared 
That © values vary with Ap p, appears to be un- 


factors or quantities are functions of Ap p,. 


consider the dependence of 0 on Ap p, 


mistakable. 

It is interesting to compare the value of the 
criterion for the mixing between two liquids ob- 
tained in these tests with the criterion that may 
be associated with the formation of wind waves. 
The latter, as determined from the observations 
of Jeffreys, has been mentioned previously. The 
value found by Jeffreys was 0.215, and this is not 
very different from the mean value of the mixing 
criterion of the tests; that is, the value 0.178. 
When the very great differences in the conditions 
of these two phenomena are considered, the close- 


ness of the two values is quite surprising. 


V. Mean Velocity Distribution in Currents 


The facet, as shown earlier in this paper, that 
the criterion of mixing takes on two different 
values, one value when the Revnolds number of 
the current is below 450 and another value when 
above, may be explained on the grounds of velocity 
distributions. For the changes of the interface 
the velocities near the interface are the important 
ones The representation of the data, however, 
was based on the mean flow, and it will be illimi- 
hating to examine the variation of the mean ve- 
locity distributions in the cross section of a square 
channel for laminar and turbulent flows 

For this purpose we take a square of side 2/ 
draw an inside square of side 2 (/—s), with the 
centers coinciding, and consider the space between 


the two (fie. 7) 


s. Denoting the mean velocity in the complet: 


This space is a freme of width 


large square by uw, and in the frame only by w 
we wish to find the variation of wu with s/. HI 
the flow is turbulent, using the Blasius law of 
velocity, and supposing that the wall shear i 
constant along the entire perimeter, it is a simpli 


matter to show that 
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Values of this relation are shown in figure 7. 
Accordingly, when the flow is turbulent, the 
velocities near the interface are not much dif- 
ferent from the mean velocity of the entire section. 
In the actual conditions of the tests, the differ- 
ences are somewhat smaller than is indicated by 
the curve, since the liquid at the interface is in 
motion, and this was ignored in the derivation 
viven above 

The computation of a similar result for the con- 
dition of laminar flow could have been made from 
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0 aminar and turbulent reqgumes, 


the formula of velocity distribution in a rectangu- 


lar channel that Bateman gives [6]. However, 
computations from this particular formula are 
somewhat cumbersome, and it was decided to 
derive a simpler formula involving algebraic poly- 
nomials only 
Taking the axis of the channel as the 2 axis, 
pressures decreasing with z, the equation of motion 
Is 
op o'u , Ow . 
az * (on tax) ae 
Where p is the pressure, uniform across a normal 
section. Let 2/ be the length of the square side. 
Let u,, be the maximum velocity in the section; 
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that is, the velocity at r=0, y=0. Writing the 
variables 





u ‘) 
X, 
u“ 
& 
/ ee .. (48) 
Y 
n. 
/ J 


and introducing the dimensionless form of the pres- 
sure term, 
1 Op 


{ 4)) 
= a ee K, i 


the equation of motion now becomes 


V-xX K ») 
where 
vy: o- \ @) 
or On 


x=0 r=], l<n<l 
x0 y= | ‘<—26) | 
i =e 70 


The velocity term may be expanded in terms of 
algebraic polynomials: 


x l aX, hx, -CXq ‘a ow A 


where ax, isa particular solution of eq 49, and the 
remaining polynomials x4, xs, - are the solution 


of 


The quantities a, 6, ¢ are numerical constants 


Since the cross section is a symmetrical one, the 
polynomials x», x4. X . . . are functions of © and 
7. Also, in particular, the section is a square, 
and thus the value of the polynomials will not be 
altered when ¢ is replaced by n and 7 1s replaced 
by ¢. Subject to these conditions 


X= oT 9 
GH + — Hey 34 
Xs = f° +- n° — 28fn* — 28" n? + TOE... J 


As ax, is solution of eq 50: 


a «/4. aD 





es . oS. aka are 5 

















For the 
determination of the constants a, 6. and ¢ we 


This determines « from the constant a 


write from eq 51 


x0, ec] ¢=U 
x=0 n=l ¢ I AL 
ox /0e=—0 y= I c=0 
Accordingly, 
a+-b+e | 
2a — 4b+- 1c I 


a Ob INC 0) 


Solving these simultaneously, 


a 0.8456, 
h O.L5DS1. (5S 
é OOO30S 


The relation between the pressure fall and the 
maximum velocity now is 


Op . ul » 


and between the mean velocity and maximum 
velocity, 
u=—0.4797u 60 


j 


The values of u,'“ as determined in this analysis 
are shown likewiseiin figure 7 by the curve labeled 
“laminar.”’ Aee ordngly, in laminar flows the 
velocities near the interface are appreciably lower 
than the mean velocity of the entire section, and 
this faet will lead to a smaller value of the mixing 


criterion if it is computed on the basis of mean flow. 


VI. Nature of Assumptions in the Analysis 
for the Criterion of Mixing 


It is necessary to call attention to the aupproxi- 
mate nature of the assumptions that were made in 
the derivation of the criterion 

The first assumption is that the velocity of 
advance of the internal waves can be obtained by 
the ordinary well-known wave formula. The 
formula refers to the case of discrete densities 
and a surface of velocity discontinuity at the inter- 


face. Actually, although the condition of discrete 
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velocities on the two sides of the interface is nearly 
satisfied, on the contrary the velocities are not 
discontinuous. In fact, at the instant just pre- 
ceding the formation of nondeformable waves, the 
flow at the interface and on the two sides of it is 
viscous, and the regime is laminar. According 
to an analysis relative to the boundary layer of 
the interface when the density difference and. the 
viscosity difference of the two liquids are small, 
the velocity at the interface is nearly half that of 
the current, the lower liquid being a pool {7]. 
During some of the runs small wax beads loaded 
with lead particles were introduced. These would 
fall through the lighter liquid, but in reaching the 
lower liquid the beads would float at the interface 
The motion of the beads, whenever this was ob- 
served, showed that the interfacial velocity was 
nearly half that of the current 
19, corresponding to the current velocities 


For example, in 
run EF 
2.29, 3.53, and 4.56 em sec, the interface velocities 
noted were 1.37, 1.78, and 2.12 em sec, respec- 
tively. Thus, in a more precise derivation of the 
velocities of the interface waves, attention might 
be given to the laminar regime of the interface 
The second assumption consists of a statement 
as to the wavelengths attained at the instant the 
mixing is about to commence. The wavelengths 
for the current velocity are Maximum values con- 


sistent with the condition of stability; that is, 


rl’ 

qAp pr - 
Whenever records of wavelengths were kept, the 
recorded values appeared to be in agreement with 
the formula, as can be seen from the data in figure 
8. The wavelengths were estimated rather than 
accurately measured. This fact and also the fact 
that the 
weaken the value of the evidence shown. 


measurements are few appreciably 
The third assumption assigns to the waves a 
value of velocity propagation in conformity with 


the wavelengths; that is, 


~ 


w= U,/2. (62 


Although in many cases this relation held nearly 
true, there were noted instances when the wave 
velocity was larger, and in the vicinity of w= 0.80. 

The final assumption is in regard to the supply 
of energy to the internal waves. The source of 


supply is thought to be due to the eddies in the 
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Ae — cm 
sion of observed with computed wave- 


lengths at initial mixing 


troughs or depressions of the upper current. The 
supply is of such magnitude as to balance that part 
of the dissipation of the internal waves taking 
place in the liquid of the lower pool. Probably 
this idea is correct. However, our attempts to 
verify it by observations on the mechanism of flow 
did not prove successful. 

Whatever be the details of a new analysis, car- 
ried out with greater rigor and with proper assump- 
tions, the parametric form of the criterion of mix- 
ing will not be modified. Whatever changes are 
obtained will be in the nature of replacing the nu- 
merical constants appearing in eq 40, either by 
new constants or by a new function involving the 
densities and viscosities of the two liquids appear- 


Ing asa correction 


VII. The Amount of Mixing With Increasing 
Current Velocity 


The second purpose of the investigation was to 
determine the amount of mixing when the velocities 
of the current are above the eritical value. The 
amount of mixing, g, may be defined as the vol- 
ume of the heavier liquid that crosses unit area of 
the interface in unit time. The dimensions of the 
quantity gq are that of velocity, and the symbol 
(,, will be used for this quantity. The tests of 
the amount of mixing were made with the large 
flume 
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The experimental determinations of g, or of U",, 
are based on the following simple concepts. Let 
\° be the volume of the upper liquid. In these 
tests, Vis the volume of the lighter liquid con- 
tained initially in the reservoir below and prior to 
its introduction into the flume. The heavy liquid 
of the pool is highly colored with a darkish-blue 
water-soluble dye. Let \N be the number of 
coloring particles in unit volume of the heavier 
liquid. Let A, be the total area of the interface 
As a result of mixing the upper liquid becomes 
tinted. Let dN be the number of coloring particles 
entering into unit volume of the upper liquid dur- 
ing the time interval Af when the current velocity 
is constant for this interval. 

Because of continuity of mass 


VdN = AgqNM. 6:3 


Introducing the symbol 


AC=dN/N, (64) 
Vac , 
q A Af ~~? 


We may refer to SC as the increase in the frae- 
tional concentration of the upper liquid. 

During the conduct of tests the velocity of the 
upper current was changed to the values 1), UW, 


Us, . at the times tf, f,. ft. . Thus, 
during the time interval t;—f), the velocity of the 
current was constant and equal to u,. In general, 
during the time interval ¢,—¢,.,, the current 
velocity was u,. At the times fo, f, ft, . . ythat 


is, at the times that new velocities were established, 
samples 8, 8), 8, . . 8, Of the liquid in the 
reservoir below were taken. Let the partial con- 
centrations of these samples be Cy, C), C2, 

(’,, respectively. Thus to obtain the rate of mixing 
g», When the velocity of the current is u,, we 
consider the equation 


66 


The determination of the partial concentrations 
(’ of the samples x was effected by the use of 
photoelectric cells. Standard solutions of known 
partial concentrations of 0, 2, 4, 6, 9, and 10 
percent, using the heavier liquid from the pool as 


base, were prepared. These were introduced one 
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ata time into a rectangular glass jar with plane 
polished surfaces. The jar was interposed be- 
tween a source of light of constant illumination- 
mtensity and a photoelectric cell, and the milli- 
ampere reading of the cell excitation was noted. 
The milliampere reading was taken as a measure 
of the partial concentration of the sample in the 
jar. These measurements of the standard solu- 
tions of known partial concentrations established 
a calibration curve. Next, the milliampere read- 
ings of the samples so, s, x, Were obtained. 
By reference to the calibration curve, the partial 
concentrations (5, 4, (’, were then deter- 
mined 

The mixing-rate values as determined by the 
above process for run E-14 are given in figure 9 
The circles correspond to the photoelectric neas- 
urements. The square marks the velocity / 
which by visual observation was associated with 
the condition of slight) mixing. The triangle 
marks what we have judged as the critical velocity, 
Y The distribution of the observation points 
marked by circles can be best approximated by a 
curved line. This same condition was observed 
in practically all of the runs of the E_ series 
Nevertheless, for the purpose of effecting sim- 
plicitvy in the representation, a straight line is 
drawn, starting from the point of slight mixing, 
Y This means that, approximately, the amount 
of mixing is proportional to the increase of the 
current velocity above the velocity [ of slight 
mixing. Accordingly, 

t/t / ( 67 


The value of ©, 


from all the runs of the series E tests, made with 


the constant of proportionality, 
the large flume, are shown in figure 10. It would 
seem that the magnitude of the constant decreases 
when the relative density of the pool liquid is 
increased We shall 
effect of density 


ignore, provisionally, the 
Furthermore, there is a corre- 
lation between the velocities (7) and / such 
that, / 1.157 On this basis, 


form of the law of mixing rate will be 


the simplified 


/ 3.5 107° (1 1.157 6S 


It was unfortunate that, due to the limitations 
of the apparatus and the capacity of the pump, it 


was not possible to have current velocities higher 
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current of | 
than those emploved in the tests. It is important 
to have an idea of the limiting velocity where the 
law ceases to be linear, or nearly linear, and of 
the marked changes that may occur at the higher 
Another aspect of the 


velocities undesirable 
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ties of the two liquid 


apparatus was the shortness of the flumes. It 
Was not always possible to divert correctly the 
dragging of heavier liquid at the exit end. In 
some cases the diversions were somewhat excessive 
and in others somewhat deficient, and, therefore 
beeause of the shortness of the flume, errors wer 
made in estimating the mixing rates, (7 

It will be interesting at this poiat to inquire 
about the increased frictional force between thi 
two liquids due to the mixing. The frictiona 
force caused by the maintenance of the waves at 


the interface may be expressed as 


ldk 
U dt’ 
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where (7 is the current velocity and dE dt the 


dissipation in the waves per unit area of interface. 
Writing for the shear, 


T) fipl” ya (70) 
where f, is the coefficient of friction, it is seen that 


» ( y 
p~ 2, HF. 


ol dt (71) 


As an estimate it may be supposed that the dis- 
sipation will take on the same value as for sur- 
face waves in contact with air. The expression is 


dk —_ » 
dt 16m pw*a* d°: (2) 


Hence, 


a | 
XK 


y w\sja\’ 
f,—32"' .( -)( ): 
TAUAUIAY 
Assume that w= (7 2,a condition that was observed 
in some tests, and 
f,=S829' (wv AL) (a n)?. 74) 


This would indicate that, a@ and ¥ remaining the 
Actually, 
however, a increases with (7° Let it be assumed 


sume, f; decreases with increasing U’. 


for the present that a? is proportional to (7 and 
that A is unaltered. Then the value of the 
coefficient can be obtained by noting the condi- 
tions at the onset of mixing, and 


where a and \ are the wave dimensions at the con- 
dition of slight) mixing. Using the observed 
values A=2.0 em, a=0.2 em, v=—0.01 em?/see., 
and (°.—10 em see., we obtain 


f, — 0.0010. (76) 


Denoting next the additional shearing resist- 
anee due to mixing by 7», 


fo=pU,U, (77) 


since the liquid crossing the interface is initially 
at rest, and after crossing the interface it: takes 
on the velocity of the current. Writing rz, in the 
form 

t2=f.pU?/2, 7S) 
then 

f=20/,/U. (79) 
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Introducing the expression for U,, from eq 6S 


f.=7.0X10-*(1—1.15U,./U). (SO) 
Thus, if the mixing law as found above is valid 
also for large values of CU’, the coefficient of friction 
due to mixing is 


ty —0.0007. (81) 


These results show that for moderately small 
channels and moderately small velocities the 
friction coefficient at the interface is smaller than 
the coefficient: ordinarily obtained for smooth- 
walled channels. Perhaps this explains our failure 
to obtain reliable readings of pressure fall along 
the flumes during the tests. 


VII. Mixing in a Large Body of Water 


The bearing of the results of the present in- 
vestigation on the interfacial mixing occurring in 
large bodies of water in stratified flow must be 
discussed. 

In the laboratory experiments, where relatively 
short flumes are used and one of the liquids, the 
lower or the upper, is at rest, and the other is 
in motion, the state of the interface is one of 
discontinuity of density. In the case when the 
lower heavy liquid moves, it will be supposed that 
the flow has continued for a long time and that 
the characteristic wave front is absent. Under 
these conditions, the interfacial stability, the 
critical velocity of mixing, and the mixing for 
velocities above the critical may be studied. All 
these, however, refer to a reach that is to be looked 
upon as a type of initial length. 

In the phenomenon transpiring in natural en- 
vironments and thus involving large bodies of 
water in stratified flow, it may be assumed that 
conditions arise so that initial reaches are estab- 
lished. It is expected that the changes taking 
place in the initial length will be similar to those 
observed in a laboratory. What is not known 
definitely in this respect is the direct applicability 
of the results to be obtained in a laboratory 
investigation to the prototype magnitudes. Cer- 
tainly, however, a qualitative similarity, at least, 
must exist. 

If that is granted, the application of the labora- 
tory results must be restricted to a very short 
reach, which will be viewed as the initial reach. 
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Bevond this in the remaining reach, which will be 
of considerable length, the conditions for mixing 
and the manner of mixing will be of a different 
type, obeying different laws. For the mixing in 
the initial reach will establish in the following ad- 
joining reach a transition laver between the 
liquids, and in this transition laver the density 
will vary continuously. As an illustration of how 
this can be brought about, we may visualize the 
following situation. A current of fresh water of 
depth //7 is flowing with a uniform velocity over a 
pool of heavier liquid. In places in the initial 
length, portions of liquid coming from the lower 
pool will spread themselves in the upper current 
only gradually Let it be supposed that the 
spreading is proportional to time, this being meas- 
ured from the instant of departure from the lower 
liquid. The density gradient established for this 
case can be obtained readily. Taking the instant 
when the spreading is completed and the area 
covered is a square of sides of length /7, the con- 


eentration along a vertical Tay be represented by 
Lee 


for finite distances away from the interface 

Now, the actual law of spreading for a given 
avetual case may not be as simple as in the above 
illustration. As long as it is assumed that spread- 
ing of liquids ejected from below into the upper 
current and in the initial length is gradual, a 
qualitatively similar law for the concentrations as 
the one mentioned above will be expected But 
these distributions imply the existence of transi- 


tion lavers. Whereas the mixing in the case of 
sharp interfaces is brought about by the ejection 
of eddies at the crests of the internal waves, the 
mixing through the transition lavers must be asso 
ciated with the momentum exchange of turbulent 
motion. This is a matter to be approached using 


Richards 


therefore, a subject 


the basic ideas of the Prandtl and 
criterion for mixing and _ ts, 


outside the scope of the present investigation 


The author expresses his appreciation to K. 
Hilding Beij, Assistant Chief of the National Hy 
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I The Remainder in Linear Methods of Approximation 
sl 
By W. E. Milne ' 
The following discussion grows out of an effort to formulate a systematic treatment 
ic of the error for such diverse processes as interpolation, numerical differentiation, numerical 
2, integration, harmonic analysis, approximation by Least Squares, approximation by equating 
a moments, and other allied operations. 
First of all, a systematic and uniform procedure is exhibited by which a desired method 
of approximation may be explicitly constructed. Second, by suitable transformation the 
remainder is put in a form in which it is possible in many cases to estimate its approximate 
magnitude. And third, the theory is applied to a variety of concrete examples, and bounds 
by are obtained for the magnitude of the error 
nd I. Fundamental! Problem of Approximation is the problem that we shall call the fundamental 
I In the various types of problems to which this ae peach ong ageremw 
oi Jee - In order that the subsequent analyses can be 
: investigation applies, we are dealing with a real performed, it will be assumed henceforth that f(x 
” funetion f(r) in an interval aSvsbh, and we are has a continuous derivative of order n+-1 in a fixed 
supplied with a more or less adequate description, interval aSvib. The values of n to be selected 
vic or definition, of this funetion. For example, we will appear as we proceed. 
1 may have a table of values of f(r) furnished by 
= experimental observations, or a table of values Il. Basic Functions 
‘aad of f(r) computed from a series or an asymptotic 
expansion, or we may have the values, not of f(r) In all the types of approximation here considered 
ay itself, but of certain of its derivatives, f’(r), or the first essential idea is to employ a basic set of 
i’(ry ete. Insuch instances f(r) is not completely linearly independent known functions —u)(4 ; 
defined by the information provided. On = the mir), ..., U(r), and to use an appropriately | 
other hand, f(7) may be a particular solution of a chosen linear combination of these functions as 
given differential equation, for which no solution an approximation to the given function. So far | 
in elementary form is known. Or f(r) may be as theory is concerned any other set f(r), 7(2). ' 
the solution of an integral equation or of some . ie eee equivalent to the set {u,(r) by linear ' 
other type of functional equation. In such cases combination, might equally well be used. In ; 
the funetion f(r) may be completely, though practical numerical analysis however, it is impor- ; 
implicitly, defined by the given equation tant to select the set best suited to the numerical ' 
The essential point is that we may not have a process to be employed. . 
simple mathematical expression for fir), and vet In numerical methods the set of basic functions ; 
ve wish to obtain particular values of f(r), or more often employed than all others put together ' 
to find its integral, or to find its derivatives—in is the set of polynomials, 1, 4, a, .. ., 2 : 
short to be able to use f(r just vs freely as if it Actually these appear more frequently in Various ' 
were expressed in terms of known elementary linear combination such as the binomial coefficient : 
‘unctions. ; functions ( 4 ). ( - ): . ( - ): the Legendre poly - 
In order to do these things, a common practice 0) I " ; 
< to replace f(z) by an approximate expression in nomials P(r), P(x), . . ., Py(2), the Lagrangian ' 
erms of suitably chosen known functions. How coefficients Ly (x), Ly" (4), . . . Ler(x), ete. 
» secure this approximate representation of f(r) Mathematics Department, Oregon State Colleg 
' 
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Next in 


nometrice functions, 1, cos wv, sin vw, .. COS mr, 


importance probably are the ftrigo- 
Sin md which may also occur Ih power form, a 


cos sv, sin sv, . cos™r, sin™r, or in products of 


terms such as sin (4—..)/2 
used ure ¢ . I1 J), 


is an Hermitian polynomial, e77L,(2) 


Other sets sometimes 
where /7,(a 
where Z J 


It is assumed henceforth that the basic functions 


is a Laguerre polynomial, ete 


u(r), mr), . . u(r) all possess continuous 
derivatives of order »—-1 in the closed interval 


asavch, and that the Wronskian determinant 


does not vanish anvwhere in the interval a<r<h 


III. Equivalence Operators 


The basic set of approximating functions having 
been selected, the next step in solving the funda- 
mental problem of approximation is to set up the 
rule by means of which the coefficients in the 
approximating linear combination of basie fune- 
tions are to be determined We shall considet 
here only objec tive, analytic methods, as opposed 


With this 


restriction, it will be found that in a great many 


to subjective and graphical methods 


cases of practical Importance, the rule is based on 
What might be called the princi ple of equivale nee, 
We shall not attempt to formulate a completely 
shall 


citing some 


general statement of this principle, but 


instead bring out its meaning by 
familiar examples in which a principle of equiva- 
lence is employed. Such examples are 
a). Equivalence of ralue at discrete points 
That is, the value of the approximating funetion 
is made equal to the value of the given function 
Mr) at se, ry, a 

by Equivalence of line r/ elements at discrete 
points. That is, the values of the approximating 
function and its first derivative are made to agree 
with fir) and f’(a 


at assigned points. This may 


also be regarded as a limiting ease of (a) when 


pairs of potnts move into coincidence 
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c Equivalence of we ighted definite integral 


That is, the value of 
dr. 


ure made to agree with corresponding eXpressio! 
formed for the approximating function 
That 


dl Equivalence of we ighted SUNS 


the values of 


} 


are made to agree with similar expressions for 1] 
approximating function 

Evidently this list can be greatly extended 

We shall henceforth in this paper consider on! 
eases in which the coefficients are complet 
determined by some appropriate application 
the principle of equivalence 

In order to treat all such cases by neans of 
comprehensive notation, it is convenient 
introduce a set of operators Op, ),, » Oa, equi 
in number to the set of basic functions, and call 
equivale nee operators What these operators act 
ally do in any particular method of approximatio 
is determined by the rule selected for setting 
the approximation. For the examples (a), ( 


ce), and (d) above the corresponding operate 


It is desirable to employ some form of analyt 
expression for the operators O,, sufficiently gen 
to include all types of linear operators likely 
arise, but sufficiently restricted to guarantee t! 
subsequent mathematical manipulations can 
carried out. Accordingly, it is assumed het 
forth that O,f(a 
Riemann-Stieltjes integrals in the form 


ean be expressed by means 


O.ftla > > r)dq 
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in which the qijir) are suitably chosen functions 
of bounded variation in the interval asvsh. 
By specialized choices of the q we may obtain 
each of the particular operators (a), (b). (e). and 
d) listed above. For example, if g,,(7) is a step 
function with unit step at 27, while the other 
q,(£) vanish we have 


7) 


O f(r) f'? (n)dq,Ar) =f" (a 


Again if dq J wiirjde and q 0 for J l, 
2 , ke, then 


O fla wii(ryt(ryde. 


Also if qo(2) isastep function having \/+-1 steps 
of magnitude wy at zy, wy atay. . . 6, Wig at Ly, 


and if q 0, for 7 >0, then 


O, fla > wiodtlse 


From the definition of O, as a sum of integrals, 
tis clear that 
1) The operation O, Is distributive with respect 


to addition, i. e 
OD) tila hola | Oflu 0 hola 


O,+ On) Ota Ota 


» 


é If .1 is constant 
OLA) | = AO Fla 


$3) The result of the operation 0, on any fune- 
tion f(r) no longer contains vas an independent 
variable and insofar as the variable w is concerned 


rat be treated as a constant. 


IV. The Approximation Formula 


Suppose that the funetion f(r) has been given in 
mv of the ways suggested in section T, that an 
appropriate basic set ou, has been chosen, that 
i rule for determining the coefficients based on the 
principle of equivalence has been selected (i. e 
the equivalence operators O,, On, 2. 2 2. O, are de- 

ded upon), and that what we desire is a formula 
ypressing the result of some given operation O 


hen performed on f(r), where O satisfies the 
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conditions imposed on the QO, in section IIL. Let 
the desired approximating linear combination of 


the u, be 


in Which the «Ws are coefficients to be determined. 
Then the principle of equivalence gives the n+ 1 


equations 


O f(t) =O, >5 Ayu;(r). 


which by the properties stated at the end of 


section [I] are the same as 


Oo fla DAO le . ] 


It is assumed heneeforth that the determinant 


On, (ar) Oy lr) oe: OO uela 

Ou, (a Ou, (2) - O,u,(r) 
)) 

Onu,(r) Ou, (a Ou,(r) 


does not vanish. Clearly the n+ eq | determine 
the coefficients .1, uniquely. 

Qur objective is to express the result of the 
operation O applied to f(r) by means of the same 
operation, O, applied to the approximating 
function S.ljwj(r), together with a remainder 
term, PR, denoting the error of this approximate 
representation. That is 


Of(r)—O S Ajuj(a R. 


When the values of the .1, as determined by 
eq | are substituted in the expression for R, the 
jeft-hand member may be written in determinant 
form, and we have the general expression for the 
error term 


Ola On F(a Ota < x el 
Our) Ou la Ou la X Ol ula 
| 
Oua Ou (a Ousla ~ Oouila R 
I) 
Ohula Oui ta Oulu 7 Oiuta 


Formula 2 provides the explicit means for setting 
up any linear approximation formula when the 
desired operation MY, the basic set v,0r), and the 
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rule for the coeflicients based on the principle of 


equivalence have been selected 

It is obvious from formula 2 that R vanishes 
Whenever flu Is any linear combination of the 
basic set | u,(a Also, R vanishes if the operator, 
O, is any linear combination of the operators Op, 
0. 0 


V. The Operator R 


Since we have assumed that the operator, O, is 
expressible by means of Riemann-Stieltjes inte- 
grals just as are the O,, it is clear from formula 2 
that Ro may be looked upon as the result of an 
operation of this same ceneral type performed on 
the funetion f(r) and might be ealled O’f(7). How- 
ever, We wish to retain the letter R associated with 
the idea of a remainder term, and hence shall use 


the notation 


indicating both that /# is a remainder and also 
that A is the result of a linear operator operating 
on the funetion f(a 

This operator also obeys the distributive prop- 


erty with respect to addition so that 


Example 2 L neal element interpolation 


flu flu f(r 
| 0 
r r | 
l 
} 
l j j Qu 
J J 2) ] J 
Example ) (Juadrature formula 
‘2 
fiarjda flag / 
e a 
B—a | 
re) x 
/ 
) 
é 
I) 
p a 
. / 
B’ a 
ds 
n | 
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Moreover, as we observed at the end of seeti 
cv, 
R(lu(zx)j=0, i=0, 1, ... 2 


These two properties of P are essential in the si 
sequent developments 

A few examples showing the formation of P| / 
in specific cases will be illuminating. In each ea 
D) is the cofactor of the element in the upper lef 
hand corner of the numerator determinant 
the values wy. 2), .. z, are all distinet it m 
be shown that in none of these cases does 
vanish 

Example | Polynomial interpolation 
Here the formula is 


hla flu hla hla 
l l l l 
lJ I I / 
lf 
D, ri I 
/ / j 
fla / I 
| () 
F l 
Ri f(r) ] 
F Js 
/ 2 lu 
fig 
R\f 
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Example 4 Taylor's series with remainder 


ha h(a hg 
l l (0) 
| Ir J l 
) 
l I _ Dy 
I rs, na 


Example 5 Least Square approg smation., 


Since any linearly independent set of functions 
can be made orthonormal by linear combination, 
the basie set itself may be assumed to be ortho- 
normal with respect to a given interval as rsb 
Then with the abridged notation 


. *) 


| fu | flx)ula dr. ete., 


* . 


the formula is 


Ae S tu S fu ede a Stu, 
Mola | 8) () 
0) l 0) Ri fir)] 
dD / J i J 
Mu I 0) ) ae l 


\ Fourier series with remainder is a particular 
case of this example. So also is Least Square ap- 
proximation using Legendre Polynomials. 

The procedure outlined above evidently has 
great flexibility, since theoretically at least there 
are infinitely many ways of choosing a basic set 
of funetions, infinitely many ways of selecting the 
rule for forming the coefficients based on the 
equivalence principle, and similarly for f(7) and 
the operator O 

It should be emphasized that although eq 2 has 
been called “an approximation formula’’, it is not 
necessarily in any sense a “good” approximation. 
Only when the remainder term R| f(x)] is small 
enough to be ignored without affecting the desired 





does not vanish in the interval aS avs 6, we may 
construct a homogeneous linear differential equa- 
tion 


Liu Dut Pay Dru 

P.(2)D*u+ ... r ryu=0 
of order n+ 1 with coefficients P(r) continuous in 
asr<bh, which has the n+1 funetions u,(r) as a 


set of linearly independent solutions. 
In determinant form this differential equation is 


u" ul u e-* ul 

“ u, u“ u“ 
: l 0 

»| & u ? u 

NW 

u“ u ul u 
Example 1. For the basic set of polynomials 
how... 2", the differential equation is simply 


Liu Dr 'u=—0 


Example 2. For the basic trigonometric fune- 
tions, 1, cos sz, sim wv, . . cos mr, sin ma, the 


equation is 
D( DP +1) (DP? 4-27) 2... (PF? +m?)u=0 


Example 3 For the basic hyperbolic set ol, 
cosh wv, sinh vw, .... cosh mar, sinh mer, the 


equation is 


Dy DP —1)( PP —2?) 2. . (LF —m*)u=0 
accuracy do we call it a good approximation 
VI. The Basic Differential Equation Example 4 For the basie exponential set 
oe. . _e"* the equation ts 
Since the functions w,(r7) have continuous 
derivatives of order +1 and the Wronskian, W, D—1)(D—2) ... (D—m)u=0 
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It is of interest to note that in each of these four For the functions ¢ Ita 
examples the differential equation has constant 

coefficients. In general, of course, the coefficients g(r, s é “(— sy" n 
are funetions of z It may also be noted that in 


the first three examples above the differential For the functions ¢ L , 


equation ts self-adjoint g(a, s ‘ (y x)" ln! 
The general solution of the nonhomogeneous 
differential equation For the set e7, e”, , Pe. 
L ul Ol), gq ms f I ‘ n! 
Is . 
: VIII. Expression for the Remainder 
Mu yt u i + o GL. ds ) : , 
: Let f(r) be a function with a continuous deriva 
tive of order n+ 1 in the region under consideration 
in Which the C's are arbitrary constants and and set up the nonhomogeneous differential 
equation 
ul i ul i ul / 
Liu Lifts 
uf N uf S u XN 
A particular solution of this equation ts obviously 
u fir), and hence by eq 9 
u x / x uf x ® 
flu Sy u J Lig x qd r.s)ds t 
u ‘ uf ‘ uv Re ue @ ‘ 
di x } 
Mu ‘ 7] ‘ u s . 
Suppose now that we have some particular oper 
u ‘ u . u s ator, R, where R is formed for the same basic set 
asis eq 6 Apply ing the operator, Rto eq 6b we get 
in view of distributive property of R, 
uM, S u S ul ‘ 2 
> ‘ , > 
me (a wks Rl f(r] Soa Rhu (a RI | Li fis yi.) 
The essential property ol gtr.s) ts threat Since Riu I 0.7 Oo. I. nw. the first term 
on the right drops out. Now FR applied to any 
S QL (Oat as for k=0, 1, mii. function o(r) has the form 
ou" (lata x for f " a 
: Rein |l=S3 vo (nydaq ar), k=<n 
VII. The Function g(x, s) f=. ; 
For the basie set 1. a r” we find that where the @g,(7) are of bounded variation in 


a<r<h 


Consequently 


kor the baste set 1. cos vw. sim a . COs mr, 
Sin mer, hi | | Li f(s G(rs is| 
“-J 2n)! 1 “y Is Jd 
Vir, s é sin ( > ali — ae | Jt Tits q ris)¢ at q J 


For the set 1. cosh wv. sinh v. . cosh ma, By successive differentiations of the integral an 


sink mer, use of eq (5) at each step we obtain 
. F » . 
F R a 0 
Ors 4 sinh 7 Dn ' | L fis)igtr, s L fis r, S)ds 
| ( y 4 ) da : , Je Ou 4 
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hence 


nlf Lscrate as] 


= o' 
> >» L) f(s)) or! gr, s ds dq (x). 


lt is convenient to define g(r, s) by the conditions 
(94, S for rI28 ) 
ly for x “) 


Then the functions 0G(z,s)/O2t, 7=0,1,. . .. n—1, 
ire continuous in the square asa<b, ass<h, 
md vanish identically in that half of the square 
for which x-—s Hence in the double integral 
above we may replace g by 9, replace the variable 
ipper limit « by the constant limit 6, interchange 
the order of the integrations, and put the summa- 
tion under the first integral. The result is 


rl ( Li fis))g(a, s as | [" L( f(s))G(s)ds, 


where we have set 


rp e O' 


os J, ays 14, 8)dai(z) 


R g(x, S |. 


i 


We have therefore 


R f(x)] L fis G(s ds. (S) 


r which 

G(s Rigta, Ss |. 9 
Note that RP operates on g(z, s) as a function of s, 
OLs Equation S prov ides the desired expression 


or the remainder term. It should be noted that 
q 8 is merely a mathematical identity and there- 


fore inone sense tells us nothing new. Actually, 
however, eq 8 accomplishes the important step 


separating the problem into two parts, one of 
vhich, LC f(s)), is independent of the operator RP, 
ind the other, G(s), independent of the function f(r). 


IX. The Function G(s) 


For brevity the function G(s) has been written 
is a function of the single variable s, since s is 
he variable of immediate concern in the integral 
feq 8. In reality G(s) depends, directly or in- 
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es 


directly, also on all the values of z, fixed or vari- 
able, that are implied in the performance of the 
operation R. When s is larger than the largest z 
involved in the operation R{ f(2)), the function 
G(s) vanishes identically because for all such values 
of s the function g(z,s) is zero by definition. Again 
if s is less than the least z involved in R{ f(x] 
the funetion G(s) vanishes identically, this time 
because q (2,8) g(a,s) and is therefore simply a 
linear combination of the w,(x), so that 


Rig (a,s)] 0). 


It follows that the limits a, 6, of the integral in 
eq 8, may as well be replaced by 2, +o, 
respectively. 

The behavior of G(s) in the interval where it 
does not vanish can best be illustrated by a 
number of simple examples. 

Example 1. Let the basic set be 1, 7, 2°, 2°, 
let 2, 2), To, Zs, 2 be five equally spaced values of x 
with interval /, let y,—f(r,) and let 


Rif(z)\= ys — 4s +-6y, — AI + Yo Aty, 


For this case 











q rs r s)° 3! 2 8s, 
whence 
(i(s RG (7,8)]= 59 [(ay—s) 4(7,;—s 
>. 
sine ~tl  + 


ay pe 
where the notation (r—s) Means that 


((rx—s)* if r>s, 


i oo if r<s 


We now see at once that for 


rs > (is 0) 
l 
ras I4 (r(s of Jd; S 
». 
, l 
Tr. = I (iis = | [\a'g x tu N 
+). 
1, 
r.a= rT ty S Ts | } I S J S | 
>). 
[ 
Tos ry (r(s ae (C2 S } 
>. 
x S 1 Gis ) 
oe &.,? roe PED ine 


le ee ee 





Pie” 


aT PE RE 


SPO SH | R CTS = 


POLED SROs + we 


2 Ee REPO Ee ee 


The graph of this function is shown in figure 1. 





Figure 1. 


Thus it is seen that G(s) consists of four ares of 
four different 


joined together continuously and with continuous 


third-degree polynomial curves 
first and second derivatives. 


Since we started with RE f(7)]—A'y, we have 


shown that 


A‘y f "(2 Gi(s ds. 


. 


when G(s) is the function just deseribed 
Example 2 Let the basic set be 1, x, 27, 2°. 2". 
let x, 4), #2 be three equally spaced points with 


interval A, and let 


RUf(2) |= Ye— Yo 5 (Yo +444 Yo)’ 


This is equivalent to Simpson’s Rule 


For this case 








_ ‘ 
q(r.s a—s) /4! 
whence 
i] 
¥ I NX J x 
(rs(8 ‘1 
h 
: Pa Le 8 4{(7,—s (xr; s ] 
We have then 
ro=8 (; & {) 
’ J S h I s 
alien . i! 0” 3 <3I 
. J s)* Ada S 
J oo } (r(s 1! T ; 1 
<7 (i(s 0 


The graph of this function is shown in figure 2 


Here A SN 
polynomials 


consists of ares of two fourth-degree 


jomed together with third-degree 
contact 
Recalling the original definition of RUf(x)] in 


this example we see that 
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h by’ , F 
(7? ’ + 
3 (YU: Wit Yo y 


8s) G(s)ds. 


The right-hand member of this equation is there 
fore a measure of the error of Simpson’s rule in the 
special form chosen here 





Figure 2 


Example 3. In the two preceding examples we 
have started with a known operator, 2, and havi 
derived from R the function G(s Now that wi 
see something of the character of G(s) in the cas: 
of a polynomial basic set, it is clear that we could 


reverse the process —-i. e., siart with a G(s 


satisfying suitable conditions as to degree and 
continuity and derive from it the operator R 


To illustrate this procedure we start with 


G(s) =0 if mn=s, 
(7(8 7 s ‘ i if 268 7 
Gi(s 0 if s I 


Here A 8 


is a polynomial of degree six, so that 
the basic set is 1, 7, 27, 2°, zt, 2°, 2°, the differentia 


equation is 17 w—0, and 





Ri fir)] f(s) (a,—s)*(s—a2zX9)°ds 


Now G(s Gi’ (s G’’(s Oat s=7o and at s=—2 


GQ’ (xz 6h G’’’ (a bh 
GY(a 72h? G(x 72h 
G(x, 360k GO (a 360h 
Gr (2 720 CG (x 720. 


Integration by parts gives accordingly 
RUF (x)) 


72h? f’’ r f’’ (x) ]- 


6h f’’’ (2,) +-f'°"" (a) ] 
360hA[ f’ (2,) +f (a)] 
720[ f(x.) —f (20). 
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ds. 


. there 
>in the 


sles we 
d have 
hat we 
1 Case 
> could 
1 G(s 
and 
tor R 


o that 
rential 


ro) | 


searc! 


fransposing and dividing by 720 we obtain the 
ormula with remainder term 


cs Pa 
"0 (n) —f"" (40) ]4 


Daw . 
r,) —f\to) 9 LU (11) +f" (20)] 10 


h : es, 
[f’ (7) tf’ (ry)] =) I “"(8) (2) 8)" (8 — Lo) ds. 
Go Jy 


20) 
This formula might be looked upon as the first 
few terms with remainder of a sort of two-point 
favlor expansion. It can be used in successive 
pproximations to secure a starting value in the 
tep-by-step integration of a differential equation. 
ts virtue lies in the fact that while it requires only 
he first three derivatives it is as accurate as a 
lavlor’s series using the first six derivatives. 

This example illustrates that fact that G(s) is a 
ind of generating function from which the parti- 
ular operator FR associated with Gis) can be 
obtained by carrying out the integration 

{"*) (2) Gis\ds. 
lt is this property which suggested the notation 
fri 

Example 4. In certain instances itt may be 
easier to determine the contact properties by 
vhich the ares composing the function G(s) are 


cr 
Ie 


ined and to construct G(s) from these data 


rather than to compute R\(r4—s "| directly. An 
instance of this is the following: Let the basic 
wt be l,a,... , 2", and let RU f(x] represent the 


difference between f(r) and the best n-th degree 
polynomial approximation by Least Squares for 
he interval (—1, 1). We shall suppose that z 
sin this interval. A brief examination will show 
hat G(s) is of degree n in s, of degree 2n +1 in s, 


n+l 


vanishes at s=1 to the same order as (1—s) 


n+! 


anishes at s 1 to the same order as (1-+-s) 
nd if Gy and G, are the polynomials to the right 
nd left respectively of s=v, then 

ie 2 


sg)" 
(i, n! + (rp. 
2 


hese data determine Gp and G, completely and 


ter considerable algebraic reduction it is found 
at 
|] —g** i a og A ) 
pan bore etten eI atey fb 
for r<s<l. 
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The foregoing examples furnish an idea of the 


character of (Gs) when the basic set consists of 
polynomials. Similarly when the basic set is 
trigonometric we would expect to find that G(s) is 
composed of trigonometric arcs welded together; 
for the exponential case, G(s) would be made up of 


exponential ares, ete. 
X. Computation of the Error 


As has been shown, the remainder ts of the form 


Ri f(x] L(f(s)) @(s)ds. 


If a<s<b is the interval in which G(s) is not 
identically zero, we have 


R{f(x)])- MK (b—a), (10) 
in which 


M maxX Lif S))| in asc hb. 


K=max| G(s)|in axs=b. 


The inequality 10 provides a bound for the mag- 
nitude of the error in all cases. 

In many particular problems this result can be 
considerably improved. For example, if G(s) does 
not change sign in the interval (4,6) we may apply 
the law of the mean for integrals and write 

R | f(2)]=LUf(—)) G(s)ds, (11) 
Where ac&<b. Note that this is the case for 
examples 1, 2, and 3 of section IX. Since in 
those examples we have the explicit expressions 
for the several ares of which G(s) is composed, 
there is no difficulty in performing the integration 
over each are and thus we can obtain the actual 
value of the integral 


G(s)ds, 
However we now show that the value of the above 
integral may be obtained from the operator R 
without making any use of the explicit form of 
G(s). Let w(x) be any solution of the differential 
equation 


Liw)=1, 


and replace f(x) in eq 8 by w(2). The result is 


R\lw(a | | Gis)ds. 





<@ .  6fer 8 et ee 


a 


Pn 


ie ee ik ee eee SE 


ee Le 


< 
> 
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By using this result in eq 1 
RU f(mJ=—LGE)R|w(a), 
provided G(s) does not change sign 
XI. Polynomial Approximation 


For the polynomial case where 


Hence we have the theorem 


if it the basi set is l, ie . = r". and 


b the function G 8 does not chanae sign, 


ther the error is 
Rif J a Rix’ |, 


where a< Eb 


Chis result contains many well-known expressions 


for the error as special cases For example, in 


the case of polynomial interpolation eq 


duces to 


In the type of interpolation where both f(7) and 
f’(x) are given by the approximating polynomial 


at the given points, eq 13 becomes 


r } . 
Ri f(x)] ; Toy) J I ry , . 


In the case of Newton-Cotes quadrature formulas, 


eq 13 also gives the same results as those obtained 


by the use of Bernoulli polynomials. 
Example 1. In the case of Newton’s 
eighth’s rule” in the form 


Sh 
Ys Ys +3y2+3y,+yo) =RhY), 


the basic set is 1, 7, 2°, 2’, 2, and the points are 


Zo, 2 h, a Qh. a 3h 


There is ho loss inh generality if we take J 


since the formula is independent of any particular 


location on the z-axis Then 
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Pix] 3 (34 3 X2*+3 X 1*+0 ') 9h°/2 


and the error turns out to be m 


yO shy” 


Ply| “Bl » RC) Kl 
us 
Here, and elsewhere below, the argument of the pr 
derivative in the remainder has been suppressed th 
The theorem also applies to a wide variety of I 
approximate formulas that are not commonly Ct 
treated in texts on interpolation and numerical! da 
Integration and to w hich the customary methods ol 
for obtaining the error do not directly apply 
\ few such formulas, together with error terms 
are listed below. 
h h h 7 | 
l) y Y] 9 y / 19 y y) 79) 
: h h 
é ; ; } } } ij 
J J 9 (Y / iy / 
hi h yf 
pay ‘4 / LOOS00 
h _ , hey hy 
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These formulas have applications in the numet 
cal solution of differential equations 
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XII. Future Investigations 


So far, the possibilities of securing usable re- 
mainder terms by the methods of this paper have 
not been fully explored. It is, for example, un- 
known whether these methods will give practically 
useful forms of the remainder for the case of 
polynomial approximation by Least Squares, for 
more general 


the case of Fourier series, and 


Fourier type expansions. These cases have, of 


course, been extensively studied by Dunham 


Jackson and others by using a very different line 
of attack. 
Another matter deserving attention is the dis- 
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covery of general theorems by which it can be 
determined for whole classes of operators whether 
or not the associated function G(s) changes sign 
For we have seen that when we know that G(s 
does not change sign the evaluation of the magni- 
tude of the error is considerably simplified. On 
the other hand, if it is necessary to compute G(s 
in order to determine whether or not G(s) changes 
sign, the benefit of this simplification is consider- 
ably reduced. Except in a few special cases such 


as polynomial interpolation, no such general 


theorems are known, 
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